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I.  INTRODUCTION 


This  research  is  a continuation  of  a long-range  effort 
being  made  in  this  laboratory26 ’ 2? » 26 » to  provide 
a sound  theoretical  explanation  for  all  of  the  processes 
which  occur  in  flame  spectrometry  from  the  time  a sample 
solution  is  aspirated  into  the  flame  until  a signal  is  read 
out.  Because  the  present  work  is  limited  to  atomic  emis- 
sion flame  spectrometry  the  following  brief  review  will 
deal  only  with  that  phase  of  flame  spectrometry.  Vine- 
fordner  has  studied  the  relation  between  spectrometer 
slit  width  and  measured  intensity;  Vinefordner,  Mansfield, 
and  Vickers  ’ have  measured  the  temperatures  in  various 
analytical  flames  and  the  efficiency  of  atomization  of  the 

sample  solution  into  analytical  flames;  Vinefordner  and 
30 

Vickers  gave  a method  of  calculating  the  limit  of  de- 
tectability in  atomic  emission  flame  spectrometry; 
Vinefordner  and  Veillon^1  studied  the  effect  of  electro- 
meter noise  on  the  limit  of  detectability;  and  Vinefordner, 
Vickers,  and  Remington"^  described  a way  to  calculate  that 
particular  concentration  of  the  atomic  metal  species  in  the 
flame  gases  at  which  self- absorption  of  the  emitted  light 
is  sufficient  to  change  the  slope  of  the  analytical  working 
curve.  Of  course  there  are  other  workers  than  those  of 


1 


2 


this  laboratory  who  are  interested  in  this  area  and  have 

published  considerable  material  relating  directly  to  it. 

18  2 3 20 

Alkemade,  ’ Broderson,  * ^ Mavrodineanu  and  Boiteux, 

7 6 1-5  15 

Gilbert,'  Gaydon,  Sugden,  ’ have  written  articles  of 
theoretical  value. 

The  major  contributions  of  this  dissertation  to  a more 

complete  understanding  of  the  choice  of  optimum  conditions 

for  flame  emission  spectrometry  are:  (a)  those  instrumental, 

flame,  and  solution  conditions  which  produce  the  largest 

values  of  the  signal-to-noise  ratio  also  represent  the 

optimum  working  conditions;  (b)  the  signal-to-noise  ratio 

50 

expressions  presented  by  Winefordner  and  Vickers^  have 

51 

been  altered  to  consider  the  electrometer  noise"^  and  the 

8 

minimum  resolving  power  slit  width  and  to  consider  higher 
concentrations  of  sample  solution  than  exist  at  the  limit 
of  detection;  (c)  experimental  measurements  of  the  values 
of  various  instrumental  parameters  used  in  the  theoretical 
calculation  of  the  signal-to-noise  ratio  have  been  made; 

(d)  the  optimum  monochromator  slit  width  for  various  experi- 
mental conditions  was  obtained  by  calculating  and  experi- 
mentally measuring  the  effect  of  slit  width  on  the  signal- 
to-noise  ratio.  This  was  done  by  calculating  the  optimum 
slit  width  from  a derived  equation  and  also  by  the  graphical 
plots  of  calculated  and  experimentally  measured  signal-to- 
noise  ratios  against  slit  width.  This  plot  showed,  not 
only  the  size  of  the  optimum  slit  width,  but  also  the 


3 


sensitivity  of  the  measurement  to  variation  in  slit  width; 

(e)  verification  of  the  derived  expressions  used  by  comparing 
experimentally  measured  and  theoretically  calculated  values 
of  signal-to-noise  ratio,  and  also  by  comparing  the  experi- 
mentally measured  and  theoretically  calculated  shapes  of 
the  signal-to-noise  ratio  versus  slit  width  curves;  (f) 
a study  of  the  effect  of  variations  of  instrument,  flame, 
and  solution  parameters  on  the  signal-to-noise  ratio  versus 
slit  width  curve;  (g)  because  the  theoretical  calculations 
are  tedious  and  time  consuming,  Fortran  programs  were 
prepared  to  obtain  the  calculated  values  of  signal-to- 
noise  ratios  for  various  instrument,  flame,  and  solution 
parameters. 


II.  INSTRUMENTATION 


The  apparatus  used  in  this  study  included  a Beckman 

DU  spectrophotometer  (Beckman  Instruments  Inc.,  Fullerton, 

California)  with  a flame  attachment;  a DU  power  supply  unit, 

and  an  externally  attached  readout  system  with  meter  or 

recorder;  also  a Jarrell-Ash  0.5  meter  Ebert  mounted  grating 

mounted  grating  monochromator  (Jarrell-Ash  Co.,  Newtonville, 

17 

Mass.)  using  an  ac  amplifier  system  and  a 10  mv  Sargent 
(Model  TR)  recorder  (E.  H.  Sargent  & Co.,  4-64-7  West  Foster 
Ave.,  Chicago,  111.  60650).  The  amplifier  of  the  Beckman 
DU  was  arranged  so  that  the  output  signal  could  be  fed 
into  (a)  the  regular  readout  system;  (b)  into  an  external 
ammeter,  or  (c)  into  a Heathkit  EUW-20A  Servo  Recorder 
(Heath  Co.,  Benton  Harbor,  Michigan).  The  system  contain- 
ing the  external  ammeter  also  contained  a helipot  zero- 
suppress  circuit  with  a 1000  unit  scale  and  1-1/2  volt 
battery.  The  mirror  in  the  flame  attachment  housing  was 
covered  during  the  measurements  so  that  the  monochromator 
would  not  collect  reflected  light  emitted  by  the  flame. 

(In  this  way  the  area  of  the  collecting  lens  or  mirror,  A, 
divided  by  the  square  of  the  focal  length  of  the  lens  or 
mirror,  F,  would  be  an  accurate  measure  of  the  number  of 


4- 


5 


steradians  (solid  angles)  of  light  from  the  flame  that 
enter  the  monochromator.)  The  condenser  in  parallel  with 
the  phototube  load  resistor  which  established  the  frequency 
response  band  width  value  (Af)  of  the  amplifier-readout 
system  for  the  DU  was  disconnected  for  this  study. 

The  same  1P28  photomultiplier  tube  that  was  used  on 
the  DU  was  also  used  when  studies  on  the  Jarrell-Ash  mono- 
chromator were  carried  out.  The  Jarrell-Ash  electrometer 
amplifier  high  voltage  power  supply  (Model  26-770)  was 
used  for  these  measurements.  The  electrometer-recorder 
combination  also  provided  a direct  means  of  obtaining  the 
photoanodic  current  in  amperes.  A calibrated  Eppley 
Laboratory  (Newport,  E.I.)  standard  light  source  was  used 

to  measure  absolute  intensities  of  the  flame  background. 

-2  -1  -1 

The  spectral  radiance  in  microwatts  mm  ster  mp.  of 
this  source  was  known  at  intervals  of  10  mu  over  the  wave 
length  range  of  250  mu  to  750  mu.  These  values  were  con- 
verted  to  watts  cm  ster  mju  for  the  purposes  of  this 
study. 


III.  THEORY 


The  basic  expression  for  the  photoanodic  signal  in 

amperes  in  atomic  emission  flame  spectrometry,  as  derived 

30 

by  Winefordner  and  Vickers  is 

i = TfIWH(A/F2)n  , (1) 

where  is  a term  to  account  for  the  possibility  that 

the  monochromator  may  not  be  set  exactly  at  the  wavelength 
of  the  emitted  spectral  line.  This  is  of  importance  when 
the  spectral  line  of  interest  is  a multiplet  which  cannot 
be  resolved  by  the  monochromator.  In  the  case  of  a single, 
sharp,  isolated  spectral  line,  with  the  monochromator  set 
at  the  line  center,  the  value  of  ^ is  unity.  In  this 
study,  ^ is  experimentally  made  to  be  unity.  The  parameter"^ 
is  the  photosensitivity  factor;  that  is,  the  current  in 
amperes  produced  at  the  anode  of  the  detector  for  each  watt 
of  radiant  power  striking  the  photocathode.  The  term  T^ 
is  the  transmission  factor  which  accounts  for  the  loss  of 
emitted  light,  due  to  reflection  from  the  lenses,  etc., 

within  the  monochromator,  and  I is  the  total  radiant 

-2  -1 

intensity  (integrated  intensity)  in  watts  cm  ster  of 
the  emitted  spectral  line.  The  parameters  V and  H are  the 
width  and  height,  respectively,  in  cm  of  the  monochromator 
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2 

slit  opening,  and  A/F  is  the  number  of  steradians  viewed 

by  the  spectrometric  system  as  long  as  the  effective  slit 

50 

aperture  is  filled  with  radiation.  The  parameter  n^ 

accounts  for  more  than  one  multiple  of  the  solid  angle 
2 

A/F  being  focused  on  the  slit  aperture  by  the  use  of 

suitable  entrance  optics.  In  the  experimental  work  done 

here,  the  mirror  in  the  flame  housing  of  the  Beckman  DU 

was  blocked  off,  and  so  the  value  of  n in  the  theoretical 

treatment  is  assumed  to  be  unity. 

In  order  to  use  Equation  (1),  it  is  necessary  to  make 

certain  assumptions.  A complete  listing  of  these 

assumptions  is:  (a)  /£  = 1;  (b)  n = 1;  (c)  the  entrance 

optics  are  sufficiently  well  aligned  so  as  to  assure  that 

the  radiation  entering  the  monochromator  completely  fills 

the  collimator,  and  (d)  entrance  and  exit  slits  are  the 

same  size.  For  the  experimental  system  used  in  this  study, 

these  four  assumptions  are  valid  (see  above  discussion). 

50 

The  equation  given  by  Vinefordner  and  Vickers  for 
the  total  root-mean-square  (r.m.s.)  fluctuation  (photo- 
anodic)  current,  Ai^,  is 

A%  “ (dip2  + Aic2  + Aie2  Af)1/2  » (2) 

where  Ai  is  the  r.m.s.  noise  (photoanodic)  signal  due  to 
o 

flame  background  flicker,  Ai  is  the  phototube  noise,  and 

sr 

Ai  •i/aF  is  the  electrometer  r.m.s.  (photoanodic)  noise 
0 ' 
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signal.  Because  Aic,  Aip  and  Aig  l/Af”1  are  independent  noise 

signals,  they  are  added  quadratically.  In  the  expression 

used  by  Winefordner  and  Vickers^  for  Aip,  it  is  assumed 

that  only  the  shot  effect  noise,  which  depends  on  the  photo- 

anodic  dark  current  through  the  phototube,  is  significant; 

that  is,  Johnson  noise  due  to  random  electron  motion  in 

the  phototube  load  resistor  is  negligible.  Also  for  the 

studies  in  this  dissertation  the  electrometer  noise 

Ai  Af1//2  is  found  to  be  negligible  (see  Section  IV  A. 7) 

6 

comnared  to  aI  and  Ai. . The  phototube  noise  current,  Ai  , 

p C sr 


is 

Aip  = C2ec  BM  Af (id+ic)]1/2  , (3) 

where  e is  the  charge  on  an  electron  in  coulombs,  B is  a 
c 

2 

constant  approximately  equal  to  1 + 1/G  + 1/G  (where  G 


is  the  gain  per  stage  of  the  photomultiplier  tube),  M is 
the  total  amplification  factor  of  the  photomultiplier  tube, 
Af  is  the  frequency  response  band  width  in  sec  ^ of  the 
amplifier-readout  circuit;  is  the  photoanodic  dark  cur- 
rent in  amperes,  and  i is  the  photoanodic  signal  in 
amperes  produced  by  the  intensity  of  the  flame  continuum. 
2?he  continuum  signal,  ic,  is  given  by 

ic  = YTfIc  s WH(A/F2)  , W 

where  s is  the  spectral  band  width  in  mu,  and  IQ  is  the 
watts  cm""2  ster-1  (or  steradiancy)  per  mu  of  spectral  band 
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width,  of  the  continuum.  Therefore,  I s is  the  total 
steradiancy  of  the  continuum  that  enters  the  instrument. 

The  expression  used  for  hi  is 

Ai  » 7 T~  Al.  s WH(A/F2)Af1/2  , (5) 

C X G 

where  Al  is  the  r.m.s.  fluctuation  in  the  steradiancy  (or 
c 

noise)  per  mu  of  spectral  "band  width  of  the  continuum,  and 

therefore  Al  s would  he  the  total  r.m.s.  noise  due  to  the 
c 

continuum  flicker  entering  the  instrument.  The  total  noise, 

50 

Airp  , is  therefore  given^  by 

AiT  = [2ecBMAf(id  + yTfWH(A/F2)Ics) 

+ {yTfWH(A/F2)Alcs  ]2Af]1/2  . (6) 

The  assumptions  involved  here  include  those  mentioned 
earlier  concerning  Equation  (1)  and  also  the  following 
assumptions.  The  signal  due  to  the  emitted  light  from  the 
metal  atoms  in  the  sample  near  the  limit  of  detection  is 
not  significant,  and  there  is  no  significant  noise  contri- 
bution from  either  the  electrometer  noise  or  the  thermal 
emission  flicker. 

However  the  study  is  not  limited  to  the  very  low 
sample  concentrations  near  the  limit  of  detection,  and  so  a 
noise  contribution  from  the  r.m.s.  fluctuations  of  the 
sample  thermal  emission  must  be  included,  that  is,  Ai.  This 
fluctuation  is  caused  in  the  same  way  as  the  fluctuations  in 
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the  continuum  intensity  and  therefore  cannot  he  treated  as 
a source  of  noise  that  is  independent  of  Ai.;  that  is, 
noise  signals  Ai  and  Ai  add  linearly.  Vinefordner  and 

G 

30 

Vickers  discussed  the  linear  addition  of  dependent  noises 
and  the  quadratic  addition  of  noises  which  are  independent 
of  each  other. 

The  spectral  hand  width,  s,  is  approximately  equal  to 

30 

the  sum  of  two  terms;  one,  s , which  is  equal  to  the  re- 
ciprocal linear  dispersion,  R^,  multiplied  hy  the  slit 
width  (W)  in  cm  and  also  a constant  term,  s , which  is 
equal  to  reciprocal  linear  dispersion,  R^,  multiplied  hy  a 

constant  slit  width  (W_)  in  cm.  Thus 

m 

s = s + s =*  Rj (W  + Wm)  (7) 

m c dN  m/ 

The  minimum  resolving  power  slit  width,  Wm,  represents  the 
minimum  width  to  which  the  slits  can  he  closed  and  still 
maintain  a properly  focused  image  on  the  dispersing  element. 
The  minimum  resolving  power  slit  width,  is  due  to:  (a) 
the  improper  focus  of  optics  within  the  instrument  leading 
to  an  inescapable  aberration  or  fuzziness;  (b)  diffraction 
effects  at  an  infinitely  narrow  slit  width  (i.e.,  the  slit 
at  very  narrow  widths  acts  as  a transmission  diffraction 
grating) . 

Therefore  Equation  (6),  with  the  assumptions  and  in- 
clusions mentioned  above,  is  then  given  more  exactly  hy 


AiT  = [2ecBMAf  [id+rTfV/H(A/F2)[IcRd(W+Vm)+I]} 
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+ 


2 1/2 


{yTfVH(A/F2)[AlcRd(V+Wm)+Al]  } Af]  . (8) 

Equation  (8)  can  be  simplified  by  combining  tbe  constant 
terms  and  substituting 


Tbe  terms  Al  and  AI  can  be  further  simplified.  Equation 

o 

(1)  shows  that  the  signal,  i,  is  directly  proportional  to 


duce  fluctuations  in  i,  Ai,  or  noise.  Fluctuations  in  I 
are  amplified  to  the  same  degree  that  I itself  is  amplified. 
Consequently,  if  the  fluctuations  of  I produce  a steady 
"white"  noise  superimposed  on  a steady  dc  signal,  the  ratio 
of  r.m.s.  noise-to-signal  should  remain  constant  at  various 
states  of  amplification  because  a change  in  the  amplifica- 
tion of  one  is  accompanied  by  an  equal  change  in  the 
amplification  of  the  other. 


K,  = 2e  BM 
d c 


(9) 


and 


Kc  = 7TfH(A/F2)  . 


(10) 


the  intensity  I;  consequently,  fluctuations  in  I,  AI,  pro- 


Al/I , 


(ID 


and 


AI  » ^1 

A similar*  situation  exists  between  IQ  and  AIQ  so  that 


(12) 


(13) 
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It  was  mentioned,  earlier  that  the  r.m.s.  fluctuations 
represented  by  Al  and  Al  are  not  indepdnent  of  each  other 
and  are  probably  caused  simultaneously  by  the  same  flame 
variations.  It  is  reasonable  to  assume  that  the  flame  is 
in  thermal  and  chemical  equilibrium  and  homogeneous  over 
the  region  being  viewed  and  that  these  flame  variations 
would  cause  the  same  degree  of  fluctuation  in  each  intensity 
simultaneously , and  so 

/ = jf>  . (14) 

Such  an  assumption  was  made  for  the  theoretical  treatments 
of  this  study.  When  Equations  (9,10,12,13  and  14)  are  sub- 
stituted into  Equations  (1)  and  (8),  then 


i = K VI 
c 

and 


(15) 


‘dlW'VdH1*11!  + »0W^)2 


[I0Rd(W+Wm)+I]2  ]1/2Af1/2  . (16) 

Using  Equations  (15)  and  (16),  the  signal-to-noise  ratio 
is  then  given  by 


i/AiT  - K0VI/[Kd  [id+KcV[IcRa(W+Wm)+I]j 


1 

(K  Wj^)2[I0Sd(W+UII1)+I]2  ] Af1/2  . (17) 


The  derivation  of  the  expression  for  the  optimum  slit 

30 

width  has  been  given  by  Vinefordner  and  Vickers.^  To  find 
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the  optimum  slit  width  Equation  (17)  is  differentiated  with 
respect  to  W,  and  when  the  derivative  is  set  equal  to  zero 
the  optimum  value  of  V is  found;  that  is,  W . According  to 
Vinefordner  and  Vickers, ^ the  differentiation  step  is 
easier  if  the  square  of  the  reciprocal  of  Equation  (17)  is 
differentiated  with  respect  to  V,  and  then  W is  optimized 
by  setting  the  derivative  equal  to  zero.  Performing  these 
steps,  and  substituting  VQ  for  W to  indicate  that  particular 
value  of  W found  to  be  the  optimum  value,  the  following 
expression  results 

(2afAI02Rd2K.2)V0*  + (2AfaIc2Rd2WmK02)U03  + 

(2AfAIoRd(()Kc2I)W05  (18) 

-CKdAfIcRdVmK0)W0  - (KdAfK0I)U0  - (2KdAfid)  = 0 . 

Equations  (17)  and  (18)  are  the  results  of  the  mathe- 
matical derivation  phase  of  this  study.  They  are  incorpora- 
ted into  a Fortran  program  (described  in  greater  detail  in 
Appendix  3)  designed  to  calculate  values  for  the  signal-to- 
noise  ratio  as  a function  of  integrated  intensity  values,  I, 
and  arbitrary  slit  width  values,  V,  after  evaluation  of  the 
various  instrumental  and  flame  compositional  parameters 
appearing  in  the  equations.  The  values  for  I are  obtained 
from  the  results  of  other  Fortran  programs  to  be  described 
later;  these  values  for  I are  ultimately  functions  of  the 
type  of  element  and  line,  the  sample  concentration,  the  flame 
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temperature,  and.  flame  type.  Therefore,  the  calculated 

values  of  signal-to-noise  ratio  and  of  optimum  slit  width 

can  he  related  to,  and  are  functions  of,  those  parameters. 

It  was  necessary  to  derive  suitable  expressions  for 

calculating  the  values  of  I as  a function  of  element  and 

line,  sample  concentration,  flame  temperature,  and  flame 

type.  The  general  pattern  for  such  derivations  was  avail- 

24- 

able  in  the  work  of  Vickers  in  which  he  derived  such  an 

expression  for  I for  the  element  sodium.  This  study  extends 

those  derivations  to  apply  to  the  alkaline  earths.  For  the 

case  of  the  alkaline  earths  and  the  alkali  metals,  it  is 

necessary  to  consider  the  metallic  compound  species  that 

can  be  present  in  significant  amounts  in  the  flame  gases. 

The  general  approach  to  the  derivation  of  an  expression 

for  I is  to  start  with  the  molar  concentration  of  the 

sample,  C,  and  find  the  total  partial  pressure,  P^,  in  the 

flame  gases  of  the  aspirated  metal  in  all  of  the  forms  in 

which  it  may  be  present  (see  Appendix  2 for  relationships 

between  C and  P^).  Then  the  partial  pressure,  PM,  of  only 

3 

the  atomic  metal,  M,  and  the  number  of  metal  atoms  per  cur 
of  flame  gases,  N,  are  calculated.  From  the  N value,  it  is 
possible  to  calculate  the  integrated  intensity  and  atomic 
absorption  coefficient  of  metal  M in  the  flame.  By  means 
of  Fortran  programs  (see  Appendix  2)  designed  for  this  study 
the  intensity  value,  I,  and  the  atomic  absorption  coefficient, 
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k°,  are  calculated.  The  former  is  used  in  the  calculations 
of  signal -to-noise  ratio  for  atomic  emission  flame  spectro- 
metry. The  latter  is  not  used  specifically  in  this  study 
but  will  be  useful  to  a study,  similar  to  this  one,  on 
atomic  absorption  flame  spectrometry. 

In  Appendix  1,  a discussion  of  the  mathematical  ex- 
pressions and  the  Fortran  program  for  the  calculation  of 
values  as  a function  of  flame  temperature,  flame  type, 
aspiration  rate,  and  molar  concentration  of  sample  is  given. 
The  results  of  this  program  are  printed  out  in  tabular  form 
where  flame  temperature  and  sample  concentration  are  the 
two  indexing  dimensions  of  a table  of  P^  values.  There  is 
a separate  page  (and  table)  for  each  gas  flow  ratio  of  each 
different  combination  of  fuel  type  and  oxidant  gas.  These 
PT  values  are  independent  of  the  type  of  metal  in  the  sample 
solution;  the  coverage  is  broad, and  both  extrapolation  and 
interpolation  are  easily  possible  in  the  use  of  the  tables. 
The  tables  produced  by  this  program  should  contain  informa- 
tion of  value  for  a variety  of  flame  studies. 

In  Appendix  2,  the  mechanics  and  format  of  the  Fortran 
program  used  to  calculate  values  of  the  integrated  intensity, 
I,  as  a function  of  element  and  line,  P^,,  flame  temperature, 
and  flame  type  are  given.  In  order  to  obtain  integrated 
intensities,  various  parameters  must  be  known  (see  Appendix 
2).  To  calculate  the  partial  pressure,  P^,  of  the  atomic 
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metal  species  in  the  flame  gases,  several  assumptions  must 
be  made.  It  must  be  decided  which,  of  the  possible  compound 
and  ion  species  of  the  metal,  will  be  present  in  the  flame 
gases  in  large  enough  quantities  so  as  to  subtract  signifi- 
cantly from  the  concentration  of  the  atomic  metal  species; 
that  is,  the  only  species  capable  of  emitting  light  at  the 
particular  spectral  line  frequency  being  studied.  The 
decision,  as  to  whether  or  not  a particular  species  is 
competing  seriously  with  the  metal  species  in  concern,  is 
based  upon  consideration  of  the  concentration  of  the  re- 
actants, the  size  of  the  formation  constant  at  the  flame 
temperature,  and  the  mechanics  of  the  formation  of  a mole- 
cule of  that  species  in  a gaseous  reaction.  Values  of 
dissociation  constants  over  the  range  of  flame  temperatures, 
for  a variety  of  species,  are  required  in  order  to  select 
properly  those  species  to  include  in  the  equation  for  PM* 

The  JANAF  tables16  provided  some  of  the  dissociation  con- 
stants; some  are  calculated  by  means  of  an  equation  given 

20 

in  Mavrodineanu  and  Boiteux;  and  some  of  them  are  taken 
from  the  tables  by  Vickers.2^  Tables  1 through  5 list  those 
values  of  the  dissociation  constants  of  the  species  assumed 
to  be  present  in  the  flame  for  each  of  the  five  elements 
studied.  Information  as  to  the  possible  gaseous  reactants 
present  in  the  flame,  and  their  partial  pressures,  is  also 
required.  This  information  is  obtained  from  tables  by 
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Hottel11  and  Zaer.^  Table  6 lists  the  values  for  the 
partial  pressures  used  in  these  calculations.  In  the  case 
of  barium,  strontium  and  calcium,  the  monohydroxide  is  con- 
sidered to  be  completely  dissociated  at  flame  temperatures. 
This  is  based  on  trends  in  the  alkaline  earth  family  and  on 
statements  by  Lagerquist  and  Sugden.  ’ * Because  the 

program  is  intended  to  be  also  used  with  other  elements  in 

which  the  monohydroxide  is  of  concern,  the  monohydroxide 
term  is  left  in  the  program  and,  for  these  three  elements, 
a suitably  high  value  of  the  dissociation  constant  (100.) 
is  used  at  each  temperature.  The  total  partial  pressure  of 
the  electrons  in  the  flame  gases,  P0,  is  assumed  to  be  due 
largely  to  the  ionization  of  the  flame  gases.  Gaydon 

shows  that  the  value  for  the  electron  pressure  due  only  to 

ionization  of  flame  gases,  p',  is  approximately  constant  in 

_9 

all  of  the  flames  used  in  this  study,  and  is  1.  x 10 
atmospheres.  In  principle,  the  ionization  of  the  metal  will 
also  contribute  to  the  electron  pressure,  P0.  This  increase 
in  electron  pressure  is  accounted  for  in  the  theoretical 
treatment  given  below.  On  the  basis  of  the  considerations 
mentioned  above,  the  species  that  are  assumed  to  be  of 
interest  under  the  flame  conditions  studied  are;  the  mon- 
oxide, the  monohydride,  and  the  monohydroxide.  In  addition 
to  these,  the  ionization  constants  for  the  first  stage  of 
ionization  of  the  metal  atom  must  be  considered  at  the  flame 
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TABLE  6 


PARTIAL  PRESSURES a OF  THE  ATOMIC  SPECIES  OF  OXYGEN  AND 
HYDROGEN,  AND  THE  NEUTRAL  HYDROXIDE  SPECIES,  IN  THE 
FLAME  GASES  OF  SEVERAL  FLAMES 


Species 
Flame  type 

0 

(atm) 

H 

(atm) 

OH 

(atm) 

Stoichiometric  ^ 

hydrogen-oxygen 
2600  °K 

3.6xlO"5 

9.2xlO“5 

3.3xlO“2 

Fuel  rich  ^ 

hydrogen-oxygen 
H/0=4.,  2600  °K 

2.3x10"^ 

2.7xl0“2 

6.2xl0“5 

Fuel  rich  ^ 

hydrogen- oxygen 
H/0=6. , 2600  °K 

1.2x10“^ 

3.1xlO“2 

3. 6x10“ 5 

Stoichiometric, 
hydrogen-air0 
2200  °K 

1 . 4x10”^ 

4.9x10"^ 

3.0xl0“5 

Stoichiometric 
ace  ty 1 ene - oxyg  en 
2700  °K 

6.3xlO“5 

3.8xlO“5 

4,6xl0“2 

Fuel  rich  c 

acetylene-oxygen 
2700  °K 

3. 6x1 O'5 

9.5xlO’5 

9.1xl0“2 

aPartial  pressures  given  in  units  of  atmospheres. 

^Values  taken  from  Hottel.11 

r 33 

^Values  taken  from  Zaer. 
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temperatures  being  used.  The  dissociation  constant  of  the 

metal  chloride  is  also  included  in  the  calculation  because 

the  metal  is  in  that  form  when  it  is  aspirated  into  the 

flame.  The  partial  pressure  of  the  atomic  chlorine  species 

is  taken  to  be  equal  to  PT  for  the  alkali  metals  (2  PT  for 

the  alkaline  earths)  because  the  only  source  of  chlorine  is 

24 

the  metal  chloride  fed  into  the  flame.  Vickers  studied 
the  effect  of  excess  chlorine  in  the  flame  due  to  the 
presence  of  HC1  in  the  sample  solution  and  found  that 
excess  chlorine  had  negligible  effect  on  the  value  of  N un- 
less the  concentration  of  chloride  is  abnormally  high.  The 
dimetal  oxide  and  the  dihydroxide  species  are  excluded  from 
consideration  because  there  is  not  a large  enough  concen- 
tration of  the  mono-form  of  these  species  present  to  form 
an  appreciable  amount  of  the  di-forms  of  these  species. 

The  values  of  the  ionization  constants  used  (see  Tables  1 

through  5)  were  obtained  by  using  an  equation  given  in 

20 

Mavrodineanu  and  Boiteux. 

Vickers2^*2"*  assumed  that  losses  in  the  atomic  concen- 
tration for  sodium  due  to  the  formation  of  the  monohydride 
and  the  monoxide  are  insignificant.  The  good  agreement 
obtained  between  experimentally  measured  and  theoretically 
calculated  working  curves,  based  on  this  assumption,  indi- 
cated that  it  is  a valid  assumption.  The  dissociation 
constants  for  the  lithium  monoxide  and  the  lithium  mono- 
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hydride  indicated  that  the  same  assumption  for  lithium 
may  not  he  as  valid.  Consequently,  for  convenience,  the 
Fortran  program  for  calculating  N for  the  alkali  metals 
also  included  these  two  species  among  the  species  consid- 
ered. For  the  alkaline  earths,  the  monohydride  and  monoxide 
must  be  considered.  In  Tables  7 and  8,  values  of  N for  both 
sodium  and  lithium,  calculated  by  both  assumptions,  are 
given;  that  is,  N values  when  the  monohydride  and  monoxide 
species  are  included  and  N values  when  these  species  are 
ignored. 

The  general  method  of  deriving  an  equation  for  PM  , in 
terms  of  P^,  dissociation  and  ionization  constants,  and 
partial  pressures  of  flame  gas  products,  is  the  same  as 
that  used  by  Vickers.  The  expressions  for  the  various 
equilibrium  constants,  each  containing  P^,  are  combined 
algebraically  and  solved  for  P^  in  terms  of  parameters  whose 
values  are  known.  In  this  case,  the  final  expression  is  not 
simple  to  solve,  and  so  the  use  of  a Fortran  program  and  a 
successive  approximation  technique  is  the  most  convenient 
way  to  obtain  values  of  P^  from  the  equation  derived  here. 

For  the  alkaline  earth  metals,  the  derivation  proceeds 

as  follows:  M represents  any  of  the  alkaline  ear«h  metals; 

K,  , IL,,  K,,  are  the  dissociation  constants  of  the  metal 
1 7 d 7 ■ 

chloride,  the  metal  monoxide,  the  metal  monohydride,  and  the 
metal  monohydroxide  respectively;  Pq,  Pg,  Pq^  are  the  partial 
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TABLE  8 

VALUES  FOR  LITHIUM  IN  THREE  FLAMES  CALCULATED  WITH,  AND  WITHOUT,  CONSIDERATION 

OF  THE  MONOXIDE  AND  MONOHYDRIDE  SPECIES 
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pressures  in  the  flame  gases  of  the  neutral  species  oxygen 
(0),  hydrogen  (H)  , and  hydroxide  (OH)  respectively;  is 
the  first  stage  ionization  constant  of  the  metal;  is  the 
partial  pressure  of  electrons  due  to  ionization  of  flame 
gases  (that  is,  with  no  sample  present),  and  PQ  is  the 
partial  pressure  of  electrons  when  those  electrons  from  the 
ionization  of  the  metal  are  included  (that  is,  with  sample 
present).  It  is  assumed  that  the  partial  pressure  of 
neutral  chlorine  (Cl)  is  twice  the  value  of  P^  because  the 
aspirated  metal  chloride  is  the  only  source  of  either  metal 
or  chloride,  and  there  is  complete  dissociation  of  the 
metal  chloride.  Therefore 


Md2 

= Pm(2Pt)2/K1  , 

(19) 

PM0 

= PMPc/K2  » 

(20) 

PMH 

* pmph//k3  ’ 

(21) 

PMOH 

= pmpoh/k4  ’ 

(22) 

K5 

= PM+Pe^PM 

(25) 

However,  because  the  number  of  electrons  added  to  the 
electron  pressure  by  ionization  of  the  metal  is  simply  equal 
to  P^+  then 

• Pe  + V 


and  Equation  (23)  becomes 


29 


K5  = V(Pe  + PM+)/PM  * (25) 

By  collecting  terms  in  Equation  (25)  and  arranging  in  the 
form  of  the  quadratic  equation  where  P^+  is  the  variable, 
then 

<(ph+)2  + p;pm+  - k5pm  » o . (26) 

The  solution  of  Equation  (26)  is 

PH+  = -Pg  +C(Pg)2  + W5Pm)1/2/2  . (27) 

The  total  partial  pressure,  P^,  is  equal  to  the  sum  of  all 
of  the  partial  pressures  of  the  metal  species,  and  so 

PT  = PM  + PMC12  + PM0  + PMH  + PMOH  + PM+  * 

Substituting  Equations  (19,20,21,22,27)  into  Equation  (28) 
gives 

PT  + PM  + [PM(2PI)2/K1]  + (pmp0A2)  + (PMPHA3)  + 
(PkWV  + (-Pg  + [(Pg)2  + W5Pm]1/2/2)  , (29) 

and  clearing  Equation  (29)  of  fractions  and  rearranging 
gives 

2(K1K2K5K4  + 4 KjKjK^Pj,2  + KjKjK^Pq  + XjK^Pjj  + 

+ KlK2K5K^[(Pe)2  + 4 K5PM]1/2  ‘ 
K1K2KjK4(P(1  + 2 PT)  « 0 . (30) 

The  value  of  P^  in  this  expression  is  solved  by  successive 
approximations  using  the  Fortran  program  described  in 


30 


Appendix  2.  A similar  derivation  for  alkali  metals  differs 
only  by  the  assumption  below 

PC1  = PT  ' 

Therefore  the  expression  for  the  alkali  metals  which  corres- 
ponds to  Equation  (30)  is 

^ + K1K3K4Pq  + ^q^2^4PH  + 

K1K2K3IWPM  + KlK2KjV(pe>2  + 4 K5PM]1/2  * 

K^gKjK^CP^  + 2 PT)  = 0 . (32) 

After  values  of  P^  are  calculated,  values  of  N can  be 
calculated  from  the  P^  value  by  means  of  the  Fortran  pro- 
gram. The  values  of  N are  obtained  from  PM  simply  by  use 
of  the  ideal  gas  law  equation.  The  partial  pressure  of  the 
metal  species  is  so  low  that  the  ideal  gas  law  is  a good 
assumption,  and  so 

N = Pjyj/kT  , (33) 

where  N is  the  atoms  per  cm^,  P^  is  the  partial  pressure 
of  the  atomic  species  in  atmospheres,  k is  the  Boltzmann 
constant  in  cm^  atmospheres  °K“1,  and  T is  the  flame 
temperature  in  degrees  Kelvin. 

The  relationship  of  intensity  of  emitted  light  to 
number  of  atoms  per  cm^  is  discussed  by  Winefordner  and  co- 
workers.^0’^2  Two  expressions  are  given  relating  I and  N. 
The  symbols  in  these  expressions  are  defined  as:  h is 
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Planck's  constant  in  erg  sec;  ^ is  the  frequency  of  maxi- 
mum intensity  of  the  spectral  line  in  sec~^;  is  the 
transition  probability  in  sec-^ ; L is  the  average  diameter 
or  thickness  of  the  flame  region  being  viewed;  gu  is  the 
statistical  weight  of  the  upper  state  involved  in  the 
transition;  k is  the  Boltzmann  constant  in  units  consistent 
with  Eu  and  T;  T is  the  absolute  temperature  of  the  flame 
gases;  Eu  is  the  excitation  energy  and  Eu=hDQ  because  the 
energy  of  the  lower  state  is  zero  for  the  resonance  line;  c 
is  the  speed  of  light  in  cm  sec-'*' ; is  the  Doppler  half 

width  of  the  spectral  line  in  sec”^ ; "a"  is  the  ratio  of 
damping  broadening  to  Doppler  broadening"^  timesT^n27  no 
units,  indicating  the  broadening  of  the  spectral  line  due 
to  classical  broadening  factors  including  the  Doppler 
broadening;  B(T)  is  the  partition  function  over  all  states 

I . [HLh  «)oAtSu/10"7  4m  B(T)  ] exp(-Eu/hT)  , (34) 

which  is  valid  for  values  of  N that  are  smaller  than  the  so- 

32 

called  "intersection  point"  value  of  N,  and 

I = (hDo2/107c)[A)D  guaNAtL/ifhTB(T)]1/2 

exp(-Eu/kT)  , (35) 

which  is  valid  at  values  of  N larger  than  the  "intersection 
point"  N value.  Equations  (34)  and  (35)  are  both  required 
because  at  higher  concentrations  self- absorption  of  emitted 
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light  results  in  I being  proportional  to  the  square  root  of 

NjWhereas  at  lower  concentrations  I is  directly  proportional 

to  N.  The  "intersection  point"  value  of  N is  that  value  of 

N above  which  self  absorption  is  significant,  and  so  I is 

given  by  Equation  (35) » and  below  which  I is  given  by 

Equation  (34).  In  the  Fortran  program  in  which  I values 

are  calculated  from  N values,  a test  is  run  on  each  N value 

to  see  if  it  is  greater  or  less  than  the  "intersection 

point"  value  in  order  that  the  correct  equation  is  used  to 

calculate  the  value  of  I.  The  "intersection  point"  value  is 

calculated  by  the  equation  given  by  Winefordner,  Vickers, 
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and  Remington^ 

Ni  = 60[aA^DB(T)/AtL\o2gu]  . (36) 

The  terms  used  are  the  same  as  those  defined  for  Equations 
(34)  and  (35)  except  that  \Q  is  the  wavelength  of  the  peak 
of  the  spectral  line.  If  Ai)d  is  substituted^2  for  in  terms 
of  atomic  weight  and  temperature,  then 

= 7.15  x 10_7[167iai)oT1/2B(T)/Wa1//2  YlnTT  SuAtL\Q2]  . 

(37) 

Equation  (37)  shows  that  the  "intersection  point"  value  of 
N is  dependent  on  the  flame  temperature.  Therefore  the 
"intersection  point"  N values,  used  in  the  program  for  de- 
termining whether  the  N being  considered  is  higher  or  lower 
than  the  "intersection  point"  N,  are  calculated  by  hand  and 
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are  fed  into  the  program  as  a set  of  data  subscripted  with 

the  same  symbol  as  the  temperature  data. 

The  calculation  of  k°  values  from  N values  is  performed 

in  the  program  by  use  of  an  altered  form  of  the  equation 

29 

used  by  Vinefordner  and  Vickers,  namely, 

k°  = [2  yirflT  I2gA.  £ M 1/2/7.15xlO“7  1 -/tT87cB(T)  ]N^1/2 , 

V U.  u dl  U 

(38) 

where  | is  a factor  which  accounts  for  broadening  factors 
other  than  Doppler  broadening.  The  form  in  which  the  values 
of  U,  I,  and  k°  are  printed  out  by  the  program  is  given  in 
Appendix  2. 


IV.  EXPERIMENTAL  MEASUREMENTS 


A.  Measurement  of  Parameter  Values 

Measurements  were  made  of  the  voltages  between  the 
cathode  of  the  photomultiplier  tube  and  the  anode,  as  well 
as  the  voltage  between  the  cathode  and  ground,  for  each 
sensitivity  setting  of  the  DU  power  supply  unit  connected 
to  the  DU.  Different  sensitivity  settings  resulted  in 
different  voltages  across  the  photomultiplier  tube.  The 
circuit  diagram  in  the  instrument  manual  of  the  DU  indi- 
cated that  the  signal  from  the  photomultiplier  tube  was 
taken  off  of  the  final  dynode.  The  anode  of  the  tube  was 
maintained  at  a more  positive  potential  than  the  final 
dynode  and  was  also  more  positive  than  the  instrument 
ground  potential.  It  was  not  convenient  to  measure  the 
voltage  between  the  cathode  of  the  tube  and  the  final 
dynode  of  the  tube  because  these  measurements  were  made 
with  the  phototube  out  of  the  instrument.  It  was  assumed 
that  the  effective  voltage  applied  to  the  photomultiplier 
tube  by  the  DU  was  that  applied  between  the  cathode  of  the 
tube  and  the  instrument  ground.  Measured  voltages  between 
the  cathode  and  ground,  and  between  the  cathode  and  the 
anode  for  each  sensitivity  setting  of  the  DU  power  supply 
unit  are  given  in  Table  9. 
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TABLE  9 

MEASURED  VOLTAGE  BETWEEN  THE  CATHODE  AND  GROUND  AND  CATHODE 
AND  ANODE  OF  THE  PHOTOMULTIPLIER  ON  THE  BECKMAN  DU  SPECTRO- 
PHOTOMETER AT  SEVERAL  SENSITIVITY  SETTINGS 


Sensitivity 

Setting3- 

Potential  between 

■L 

cathode0  and  in- 
strument ground 
(volts') 

Potential  between 
cathode  and  anode 

(volts') 

Full 

640. 

720. 

9 

535. 

625. 

8 

^55. 

538. 

7 

395. 

480. 

6 

3^0. 

422. 

5 

295. 

372. 

4 

250. 

330. 

3 

210. 

290. 

2 

175. 

255. 

1 

150. 

252. 

Sensitivity  dial  settings  on  the  DU  power  supply  unit. 

^Cathode  post  of  the  RCA  IP28  photomultiplier  tubs 
socket  (the  photomultiplier  tube  was  not  in  the  instrument 
during  the  measurement;. 

cAnode  post  of  the  RCA  IP28  photomultiplier  tube  socket. 
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In  the  program  used  for  calculating  signal-to-noise 
ratio  and  optimum  slit  width  values,  the  anodic  dark  current, 
id,  the  gain  factor,  M,  and  the  amperes  per  watt  factor,  / , 
are  treated  as  constant  values.  No  allowance  was  made  for 
variation  of  their  values  to  correspond  to  the  different 
voltages  of  different  sensitivity  settings.  Most  of  the 
signal-to-noise  ratio  measurements  made  on  the  DU  were  ob- 
tained on  one  of  three  sensitivity  settings:  Full,  9,  or  8. 

Therefore,  the  voltage  applied  to  the  photomultiplier  tube 
was  taken  as  530  + 10  volts  which  was  an  average  value  of 
these  three  settings.  The  dark  current  was  measured  at  522 
volts,  and  M and  ^ were  measured  at  535  volts. 

1.  Photoanodic  dark  current  values,  i^ 

The  photoanodic  dark  current  of  the  1P28  photomultiplier 
tube  used  in  this  study  was  measured  using  the  Jarrell-Ash 
monochromator  and  the  Jarrell-Ash  dc  electrometer-readout 
system.  With  this  instrumental  arrangement,  the  signal  was 
read  directly  in  amperes.  The  signal  produced  when  a certain 
voltage  was  applied  across  the  photomultiplier  tube,  with  the 
slit  closed,  was  considered  to  be  the  photoanodic  dark  cur- 
rent flowing  through  that  tube  under  that  voltage.  In  the 
Jarrell-Ash  system,  the  signal  was  taken  from  the  anode  of 
the  photomultiplier  tube  rather  than  from  the  final  dynode 
as  in  the  Beckman  DU.  The  voltage  between  the  cathode  and 
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ground  on  the  DU  corresponded  to  the  voltage  between  the 
cathode  and  anode  of  the  photomultiplier  tube  in  the  Jarrell- 
Ash.  Measured  values  of  i^  at  various  potentials  applied 
across  the  photomultiplier  tube  are  listed  in  Table  10. 

2.  Photosensitivity  factor  values,  Y 

The  values  of  in  amperes  per  watt  were  measured  using 
the  same  1P28  photomultiplier  tube  that  was  used  in  the 
Beckman  DU  and  the  Jarrell-Ash  instrumental  system.  The 
filament  of  a standard  light  source  was  focused  by  a bi- 
convex lens  to  give  a spot  of  light  on  the  slit  opening  of 
the  monochromator  which  completely  illuminated  the  slit. 

Care  was  taken  to  see  that  the  aperture  of  the  lens  outside 
the  slit  was  less  than  the  aperture  of  the  collecting  lens 
inside  the  instrument.  Thus  the  inside  lens  was  the  one 
which  determined  the  solid  angle  of  light  that  was  eventually 
focused  on  the  grating  of  the  monochromator.  The  resulting 
signal  was  read  directly  in  amperes  on  the  recorder.  The 
slit  area  and  internal  optics  of  the  monochromator  were 
known  and,  consequently,  Equation  (1)  could  be  solved  for  Y 
which  was  the  only  remaining  unknown  parameter.  The  values 
of  Y were  measured  at  the  wavelengths  of  each  of  the  spectral 
lines  being  studied;  the  voltage  across  the  photomultiplier 
tube  was  the  same  for  each  wavelength  setting,  but  the  slit 
width  was  altered  in  each  case  to  produce  an  accurately 
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TABLE  10 

MEASURED8,  PH0T0AN0DIC  DARK  CURRENT  VALUES  AT  SEVERAL 
VOLTAGES  APPLIED^*  TO  THE  PHOTOMULTIPLIER  TUBE 


Tube  voltage 
(volts) 

Dark  current 
(amperes) 

650 

1.6  x 10"10 

609 

, c , n-io 

1.5  x 10 

522 

1.0  x 10"10 

155 

8.0  x 10-11 

aMeasured  on  the  dc  electrometer  readout 
system  of  the  Jarrell-Ash  instrumental  arrangement 
used  in  this  study. 

^Voltage  applied  between  the  anode  and 
cathode  of  the  RCA  IP28  photomultiplier  tube  by  the 
dc  power  supply  unit  of  the  Jarrell-Ash  instrumental 
arrangement  used  in  this  study. 
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measurable  signal.  Values  of  / as  measured  for  each,  wave- 
length setting  are  listed  in  Table  11. 

3.  Photomultiplier  tube  gain  values,  M 

An  attempt  was  made  to  measure  the  values  of  M,  the 
total  gain  or  magnification  factor  of  the  photomultiplier 
tube.  The  plan  was  to  short  out  the  dynodes  of  the  photo- 
multiplier tube,  to  replace  them  with  resistors  equal  to 
the  resistances  between  the  dynodes,  and  to  compare  the 
signals  obtained  with  the  dynodes  in  the  circuit  and  the 
signals  obtained  with  the  dynodes  shorted  out  of  the  circuit. 
These  efforts  to  measure  M were  not  entirely  satisfactory. 

The  signals  obtained  with  the  dynodes  shorted  out  were 
quite  insensitive  to  variations  in  slit  width  and  gave 
nearly  constant  values  at  each  wavelength  setting.  A single 
approximated  value  of  M,  very  close  to  the  value  obtained 
with  the  measurement  technique  described  above,  was  used 
in  the  program  for  calculating  the  signal-to-noise  ratios. 
This  value,  10^  (no  units)  obtained  by  measuring  the 

particular  photomultiplier  tube  used,  was  tenfold  smaller 

p ll 

than  that  used  by  Vickers  who  estimated  a value  of  M at 
10^  based  on  the  manufacturer's  claims  for  a good  photo- 
multiplier tube.  The  value  of  10^  for  M,  at  this  tube 
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voltage,  was  also  indicated  by  Lion. 

In  this  work,  it  was  assumed  that  M and  were  not 
independent  but  were  rather  directly  and  simply  related  to 


TABLE  11 


MEASURED  VALUES  OE  THE  PHOTOSENSITIVITY  FACTOR3, 
AT  SEVERAL  WAVELENGTHS 


Wavelength 

O 

(A) 

Photosensitivity 
factor  (amp.  watt-1) 

6708 

1.5 

5890 

87. 

5536 

210. 

4601 

570. 

4227 

450. 

3 A voltage  of  535  volts  was  being  applied 
the  cathode  and  anode  of  a RCA  IP28  photo- 
multiplier tube  on  the  Jarrell-Ash  instrumental 
arrangement  used  in  this  study. 
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each  other.  Therefore  this  pair  of  parameters  were  varied 
simultaneously  by  the  same  factor  as  that  used  on  a single 
independent  parameter. 

4.  Continuum  intensity  values,  I 

The  measured  values  of  the  flame  continuum  intensity, 

-2  -1  -1 

I , in  watts  cm  ster  mp  were  obtained  by  two  different 
c 

methods;  approximate  values  for  this  parameter  are  also 

7 

available  from  a graph  by  Gilbert.  The  first  method,  and 
the  method  whose  values  of  I were  eventually  used  in  the 

G 

theoretical  treatment  of  this  study,  involved  measurement 
of  the  signal  of  the  blank,  i , run  with  each  measurement 
of  a sample  signal,  i . For  each  signal  measurement,  a 

o 

base  line  (shutter  closed)  was  established,  a blank  (dis- 
tilled water  being  aspirated  into  the  flame)  was  measured, 
and  the  sample  (sample  solution  being  aspirated  into  the 
flame)  was  measured.  There  is  also  available,  in  the  data 

printed  out  by  the  programs  described  in  Appendix  1 and  2, 

-2  -1 

values  of  the  sample  intensity,  I,  in  watts  cm  ster  for 
that  particular  sample  concentration,  aspiration  rate, 
flame  condition,  and  spectral  line.  If  it  is  assumed  that 

these  theoretical  values  of  I are  valid,  then  the  ratio  of 

-2  —1 

I (intensity  of  the  sample  in  watts  cm”  ster  ) to  IQW 

-2  -1 

(intensity  of  the  flame  continuum  in  watts  cm  ster  — 
assuming  that  distilled  water  adds  nothing  to  the  intensity 
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of  the  signal  of  the  blank  but  also  serves  to  make  the 
temperature  and  shape  of  the  flame  the  same  for  the  blank 
as  it  is  for  the  sample)  is  equal  to  the  ratio  of  the  sample 
signal  to  the  water  signal.  This  proportion  is  expressed 
in  the  equation 

■ V1* 

The  values  of  the  ratio  on  the  right  side  of  Equation  (59) 
can  be  measured  on  the  recorder  chart  paper.  Using  these 
ratios,  and  the  theoretical  values  of  I,  the  absolute 
values  of  I can  be  calculated  by  Equation  (39) • The 

O W 

units  of  the  continuum  term,  I , to  be  used  in  Equation  (17) 

V 

-2  -1  -1 

are  watts  cm  ster  mp  . Values  of  I can  be  converted 

cw 

to  I„  values  by  dividing  I...  by  the  spectral  slit  width,  s. 
c cw 

The  value  of  s can  be  obtained  from  Equation  (7)  (s  = 

R,W  if  V » V ).  Values  of  R,  are  available  in  the  manu- 
d m d 

facturer's  literature  for  the  monochromator.  The  values 

obtained  by  this  method  of  measuring  I.  are  only  as  accurate 

as  the  theoretically  obtained  values  of  I. 

The  second  method  of  measuring  I consisted  of  using 

c 

a standard  light  source  on  the  Jarrell-Ash  monochromator — 
Sargent  recorder.  The  signal  in  amperes,  produced  when  the 
standard  light  source  was  focused  on  the  entrance  slit  of 
the  monochromator,  was  measured.  Also  the  signal  produced 
when  distilled  water  was  aspirated  into  a flame  placed  in 
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the  same  optical  location  of  the  standard  light  source  was 

measured.  The  ratio  of  these  two  signals  is  the  same  as  the 

ratio  of  signals  on  the  right  side  of  Equation  (39),  and.  so 

these  signal  ratio  values,  along  with  the  values  of  I taken 

from  the  calibration  table  for  the  standard  source,  make  it 

possible  to  use  Equation  (39)  to  calculate  values  of  Icw* 

The  values  of  I thus  obtained  can  also  be  converted  to 

cw 

corresponding  values  of  I by  dividing  them  by  the  spectral 

slit  width  (s)  if  W » W . In  the  case  of  the  Jarrell-Ash 

monochromator,  which  is  a grating  instrument,  a single  value 

of  is  valid  for  all  wavelengths  (as  compared  to  the  DU 

which  is  a prism  instrument).  However,  with  the  Jarrell- 

Ash  measurements,  it  was  necessary  to  use  a different  slit 

width  with  the  standard  source  than  was  used  with  the  flame 

because  of  the  large  difference  in  intensity  of  the  two. 

This  made  it  necessary  to  correct  the  flame  intensity 

reading  by  a factor  that  was  the  square  of  the  ratio  of  the 

two  slit  widths.  The  most  likely  source  of  error  in  this 

method  of  measuring  I was  in  the  mechanical  difficulties 

c \ 

of  obtaining  the  same  flame  conditions  and  optics,  for  the 

I measurement,  as  was  used  on  the  signal  measurement  with 
c 

which  the  theoretical  calculations  using  I are  supposed  to 
agree.  The  flame  used  for  the  Ic  measurement  on  the  Jarrell- 
Ash  might  not  be  in  exactly  the  same  location  of  the  entrance 
optics  as  the  standard  source  was.  The  height  of  the 
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entrance  slit  of  the  monochromator  above  the  tip  of  the 
burner  could  be  different  between  the  Jarrell-Ash,  where  I 

c 

was  measured,  and  the  DU  where  the  signal-to-noise  ratio 
was  measured.  Therefore  the  monochromator  would  view  a 
different  region  of  the  flame  which  could  have  a different 
temperature  and  different  continuum  intensity.  The  same 
burner  was  used  on  both  instruments,  but  it  was  mounted 
differently  for  each.  Therefore  some  differences  of  gas 
flow  rate  and  aspiration  rate  between  the  two  different 
mountings  of  the  burner  could  occur.  It  is  likely  that  the 
voltages  applied  to  the  photomultiplier  tube  by  the  two 
systems  differed  although  this  will  be  shown  not  to  be  a 
source  of  significant  error.  Most  of  the  errors  dealing 
with  variations  of  optics  and  flame  conditions  were  not 
present  in  the  measurements  made  on  the  DU  in  the  first 
method  because  the  same  flame,  aspiration  rate,  gas  flow 
rate,  and  monochromator  setting  that  was  used  for  the  IQ 
measurement  (distilled  water  aspirated  into  the  flame)  was 
used  only  a fraction  of  a minute  later  for  the  signal 
measurement  (sample  solution  aspirated  into  the  flame). 
However  in  the  first  method,  it  was  necessary  to  rely  on 
theoretically  calculated  intensity  values.  Consequently 
there  was  some  doubt  as  to  the  validity  of  either  method 
for  use  in  measuring  I . It  was  decided,  on  the  basis  of 
agreement  between  experimentally  and  theoretically  produced 
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curves  of  signal-to-noise  ratio  versus  slit  width,  that  the 

I values  measured  on  the  DU  gave  slightly  better  agreement, 
o 

and  so  these  were  the  ones  that  were  used  in  the  theoretical 

treatment  of  this  study.  In  Table  12  the  values  of  Ic 

measured  on  the  DU  and  on  the  Jarrell-Ash,  are  given  as  well 

n 

as  values  read  from  the  graph  of  Gilbert,  for  each  flame 
type  and  each  wavelength  studied. 

5.  Flame  flicker  parameter  values.  Ip 

The  flame  flicker  noise  was  assumed  to  be  a constant 
fraction,  }j)  , of  the  emitted  intensity.  Most  of  the  other 
sources  of  noise  (shot  effect,  Johnson  noise,  etc.)  arose 
in  the  instrument  and  were  not  dependent  on  flame  condi- 
tions. Consequently  they  did  not  increase  with  intensity 
as  the  flame  flicker  noise  did.  Therefore  at  high  sample 
concentrations,  and  high  intensity  values,  the  flame 
flicker  noise  should  be  the  predominant  noise,  and  the 
others  should  make  a negligible  contribution  to  the  total 
r.m.s.  noise,  Ai^.  The  continuum  intensity  did  not  increase 
with  the  intensity  of  the  atomic  line,  and  so  the  fluctuation 
in  the  continuum  intensity  due  to  flame  flicker  should  also 
not  increase  with  emitted  intensity  from  the  sample  atoms. 
Therefore,  at  a high  sample  intensity,  the  flicker  in  the 
flame  continuum  should  also  make  a negligible  contribution 
to  Aijp  even  though  such  noise  has  its  origin  in  the  flame 


TABLE  12 

FLAME  CONTINUUM  INTENSITY  VALUES, a OBTAINED  BY  TWO  METHODS  OF  MEASUREMENT  AND  ONE 
LITERATURE  SOURCE,  AT  FIVE  WAVELENGTHS  IN  THREE  FLAME  TYPES 
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conditions.  When  these  assumptions  are  applied  to  Equation 
(17),  those  terms  in  which  or  I appear  will  drop  out 
leaving 

i/AiT  = KcWI/Af1/2KcWI  = 1/Af 1/2 . (40) 

If  the  signal -to-noise  ratio  at  high  I values  is  known,  and 
if  Af  is  known,  Equation  (40)  can  he  solved  for  the  para- 
meter (})  . In  this  study,  four  of  the  measured  values  of 
signal-to-noise  ratio  at  different  slit  widths,  for  the 
most  concentrated  solutions  of  both  barium  and  strontium  were 
averaged  for  each  type  of  flame.  In  Table  13,  values  of 
signal-to-noise  ratio  were  averaged  for  each  flame  type  and 
calculated  values  of  If  determined  for  each  flame  are  given. 

6.  Frequency  response  band  width  values,  Af 

According  to  Equation  (17) » the  signal-to-noise  ratio 
is  inversely  proportional  to  the  square  root  of  Af.  Errors 
in  the  value  of  Af  should  cause  an  upward  or  downward  dis- 
placement of  the  signal-to-noise  ratio  versus  slit  width 
curve  but  should  not  distort  the  shape  of  the  curve.  There- 
fore, in  evaluating  the  agreement  between  experimental  and 
theoretical  signal-to-noise  ratio  versus  slit  width  curves, 
more  attention  was  given  to  agreement  in  the  shape  of  the 
curves  and  less  to  the  vertical  alignment. 

All  of  the  noises  considered  in  Equation  (16)  are 
white;  that  is,  they  are  equally  distributed  at  all 
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TABLE  13 

MEASURED  VALUES  OE  THE  SIGNAL-TO-NOISE  RATIO  AT  HIGH  CONCEN- 
TRATIONS FOR  THE  CALCULATION  OF  FLAME  FLICKER  FACTOR  VALUES 


Flame  type 

11° 

IIId 

High  concentration 

55. 

70. 

38. 

Signal-to-noise  ratio 

59. 

95. 

43. 

Values  considered 

64. 

78. 

30. 

61. 

89. 

38. 

41. 

56. 

28. 

45. 

61. 

32. 

45. 

95. 

31. 

47. 

75. 

35. 

Average  signal-to- 
noise  ratio  value 

52. 

77. 

34. 

Calculated  flame 
flicker  factor 

5.2x:10’3 

2.2xao“3 

4.8xl0"5 

aA  Af  value  of  36  cps  was  used  for  these  calculations. 
^Stoichiometric  hydrogen-oxygen  flame. 
cFuel  rich  hydrogen-oxygen  flame. 

^Stoichiometric  acetylene-oxygen  flame. 
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frequencies.  The  Beckman  DU  circuitry  includes  an  RC  filter 
network  on  the  signal  introduced  into  the  input  of  the  DU 
amplifier.  This  filter  reduces  the  value  of  Af  and  thus 
cuts  out  a considerable  amount  of  the  noise.  The  condenser 
of  this  filter  was  disconnected  for  the  signal-to-noise  ratio 
measurements  of  this  study.  The  only  filtering  action  left 
in  operation  was  the  time  constant  of  the  Heathkit  recorder. 
If  a noise  of  sufficiently  high  frequency,  such  that  the 
recorder  pen  could  not  keep  up  with  it,  should  try  to  come 
through  with  the  signal,  that  frequency  of  noise  would  not 
be  recorded.  Consequently  the  value  of  Af  was  limited  by 
the  time  it  took  the  pen  to  travel  over  that  noise  band 
distance  on  the  chart  (see  Figure  1).  The  estimated  average 
width  of  the  noise  band  at  mid-chart  was  0.28  inches.  The 
time  required  for  the  pen  to  travel  full  scale  (10  inches) 
is  one  second.  Therefore  the  value  of  Af  is  36  cps. 

This  method  does  not,  however,  account  for  the  change 
in  amplification  with  change  in  sensitivity  setting  on  the 
DU  power  supply.  If  the  input  signal  should  fall  low  on 
the  chart,  the  average  noise  signal  will  have  a width  less 
than  the  0.28  inches  which  it  would  have  had  at  the  mid- 
chart location.  Therefore  the  value  of  Af  for  measurements 
in  which  the  noise  band  width  is  less  than  0.28  inches  is 
greater  than  36  cps.  Conversely,  an  average  noise  signal 
falling  high  on  the  chart  would  be  under  a higher  amplifi- 
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cation  factor  and  would  establish  a thicker  noise  band  than 
0.28  inches  and  would  require  a longer  time  to  do  so,  and 
so  a value  lower  than  36  cps  would  result.  The  computer 
program  for  calculating  signal -to-noise  ratio  treated  the 
parameter  Af  as  a constant.  It  did  not  allow  for  variation 
of  the  value  of  Af  with  pen  position,  and  so  some  single 
value  of  Af  had  to  be  decided  upon.  The  average  mid-chart 
width  of  the  noise  band  was  chosen  to  establish  that  Af 
value  (36  cps)  which  was  used  in  the  theoretical  treatment 
of  this  study. 

To  account  for  the  variation  in  Af  with  recorder  pen 
position,  the  measured  noise  was  multiplied  by  a factor  to 
account  for  this  variation  and  to  normalize  all  noises  to 
the  middle  of  the  recorder  chart.  Then  theoretical  and 
experimental  values  of  signal-to-noise  ratios  should  be 
directly  comparable.  The  r.m.s.  noise,  as  measured,  was 
multiplied  by  a correction  factor  based  on  the  decade  of 
the  chart  paper  in  which  the  signal  carrying  the  noise 
happened  to  fall.  The  measured  r.m.s.  noise  was  decreased 
for  measurements  in  decades  lower  than  the  fifth  decade  and 
was  increased  for  signals  falling  in  decades  above  the  fifth. 
The  effect  of  the  correction  factor  was  to  give  all  of  the 
r.m.s.  noise  measurements  that  value  which  they  would  have 
had  if  they  had  all  been  measured  in  the  fifth  decade. 


51 


An  attempt  was  made  to  verify  the  assumption  that  the 
limitation  of  pen  speed  constituted  a filtering  device 
operating  in  the  instrumental  arrangement  and  was  therefore 
actually  establishing  the  frequency  response  band  width,  Af. 
A dc  power  supply  and  an  audio  frequency  generator  were  both 
connected  to  the  Heathkit  recorder.  The  dc  signal  could  be 
adjusted  so  as  to  place  the  pen  in  any  decade  of  the  chart 
paper.  An  independent  ac  noise  signal  was  imposed  on  top 
of  the  dc  signal.  Both  the  energy  and  the  frequency  of  the 
noise  signal  could  be  varied  from  zero  upwards.  With  this 
arrangement,  it  was  possible  to  determine  Af  as  a function 
of  the  fraction  of  full-scale  recorder  response.  However, 
this  instrumental  arrangement  also  revealed  that  the  pen 
speed  of  the  recorder  would  readily  pass  harmonics  of  higher 
frequencies  than  the  estimated  36  cps.  Therefore,  it  was 
impossible  to  obtain  an  absolute  value  of  Af.  As  a result, 
little  significance  was  attached  to  any  vertical  displace- 
ment of  the  theoretical  signal-to-noise  ratio  versus  slit 
width  curve  from  the  experimental  curve.  More  significance 
was  placed  on  the  curve  shape  which  was  not  affected  by 
errors  in  the  value  of  the  parameter  Af. 

7.  Electrometer  and  photomultiplier  tube  noise  values. 


With  the  shutter  on  the  DU  closed  and  the  sensitivity 
set  at  various  photomultiplier  tube  voltages,  the  total 
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noise  on  the  recorder  (pen  set  at  mid-chart)  was  due  to  the 
quadratic  sum  of  the  shot  effect  (and  other  noises  origina- 
ting in  the  photomultiplier  tube)  plus  the  electrometer 
noise  (originating  in  that  part  of  the  circuitry  following 
the  photomultiplier  tube) 

AiT  = (Aie2Af  + Aip2)1/2  , W 

where  Ai^fiiT  is  the  electrometer  noise  and  Aip  is  the  photo- 
tube noise.  An  expression  for  the  total  noise  that  includes 
the  electrometer  noise  would  be  an  alteration  of  Equation 
(16)  as  follows: 

- [ai(  Kd[id  + K0W  (l0KaCW+Vm)+I  ) ] ♦ (K0W  jy)2[I0Rd(W+Wm) 

+ I]2  + ile2j]1/2  c^2) 

If,  with  the  same  instrumental  situation,  the  voltage  across 
the  photomultiplier  tube  is  reduced  to  zero  by  turning  the 
sensitivity  setting  to  "off,"  the  recorded  noise  will  then 
be  only  the  electrometer  noise.  Then  by  adjusting  the  DU 
sensitivity  setting  at  "off,"  the  value  of  Aip  could  be 
solved  for  using  Equation  (41).  The  electrometer  noise 
signal,  AiJ&f\  was  found  to  be  negligible  compared  to  Aip 
and  so  was  eliminated  from  consideration  in  the  derivation 
of  Equation  (17). 


B.  Methods  of  Obtaining  and  Evaluating  Experimental 
Measurements  of  Signal-to-Noise  Ratio 
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Three  different  sample  concentrations  differing  by  ten- 
fold and  ranging  over  a hundredfold  were  chosen  for  each  of 
the  five  elements  studied.  However,  the  same  series  of 
concentrations  was  not  used  for  all  elements.  Each  element 
was  studied  on  the  DU  instrumental  arrangement  in  order  to 
determine  the  hundredfold  concentration  range  over  which 
signals  could  be  measured.  The  lowest  concentration  for 
most  of  the  elements  studied  was  somewhat  above  the  limit 
of  detectability.  This  instrumental  arrangement  had  a 
lower  sensitivity  due  to  the  blocked  mirror,  and  so  there 
was  difficulty  in  covering  the  full  hundredfold  concentra- 
tion range  in  some  cases.  The  concentration  ranges  and 
spectral  lines  that  were  used  are: 


Barium 

10 

ppm 

to 

1000 

ppm 

5536  A 

Calcium 

0.1 

ppm 

to 

10 

ppm 

4227  A 

Lithium 

10 

ppm 

to 

1000 

ppm 

6708  A 

Sodium 

0.1 

ppm 

to 

10 

ppm 

5890  A 

Strontium 

1.0 

ppm 

to 

100 

ppm 

4601  A 

All  of  the 

spectral 

lines  studied, 

except  sodium, 

single,  sharp,  and  isolated  as  required  by  the  assumptions 
made  in  deriving  Equation  (17).  In  the  case  of  the  sodium 


doublet,  5890-5896  A,  the  monochromator  wavelength  setting 
on  the  DU  was  maximized  with  respect  to  the  signal,  and  the 
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setting  fell  between  the  two  lines.  In  this  case,  of  course, 
the  ^ factor  described  in  Equation  (1)  does  not  have  a value 
of  unity.  However,  this  error  was  ignored  in  the  calcula- 
tions because  it  was  assumed  that  the  error  in  the  calculated 
signal -to-noise  ratio  introduced,  was  negligible  compared  to 
some  of  the  experimental  errors  introduced  into  the  measured 
signal-to-noise  ratio.  This  difficulty  did  not  arise  during 
the  use  of  the  Jarrell-Ash  monochromator  which  was  capable 
of  isolating  the  separate  lines  of  the  sodium  doublet.  In 

O 

that  oase,  the  5890  A line  was  measured  because  it  is  the 
stronger  of  the  two. 

Standard  solutions  of  all  metal  salts  as  chlorides  were 
prepared.  As  long  as  the  metal  salt  was  nearly  completely 
dissociated  in  the  flame,  it  made  little  difference  which 
salt  of  the  metal  was  used  in  preparing  the  sample. 

Each  of  the  three  concentrations  of  each  of  the  five 
elements  was  measured  at  a series  of  slit  width  values  in 

i 

each  of  three  different  flame  types:  stoichiometric 

hydrogen-oxygen  (H2/O2  = 2),  fuel  rich  hydrogen- oxygen 
(H2/02  = 4),  and  stoichiometric  acetylene-oxygen  (C2H2/O2  » 
2/5). 

The  aspiration  rate  was  measured  prior  to  each  series 
of  recorder  readings  for  each  concentration  of  each  element. 
The  aspiration  rate  was  measured  by  recording  the  time 
necessary  to  aspirate  4.00  ml  of  water.  Three  or  four  such 
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measurements  were  averaged  for  each,  reported  measurement  of 
the  aspiration  rate. 

The  signal  and  noise  measurements  were  made  on  the  DU 
under  the  following  instrumental  conditions:  the  zero- 

suppress  on  the  DU  power  supply  was  on  the  number  2 posi- 
tion; the  sensitivity  dial  on  the  DU  was  in  full  clockwise 
position;  the  selection  switch  on  the  DU  was  set -at  1;  the 
zero-suppress  on  the  external  ammeter  was  set  at  midscale 
(500);  the  damping  on  the  recorder  was  turned  off;  both  the 
zero  position  dial  and  the  calibration  dial  on  the  recorder 
was  turned  off;  both  the  zero  position  dial  and  the  cali- 
bration dial  on  the  recorder  were  in  full  counterclockwise 
position;  the  mv  range  selection  switch  on  the  recorder  was 
set  at  10.  The  above  settings  were  maintained  throughout 
all  of  the  signal  and  noise  measurements  on  the  DU.  The 
three  controls  that  were  varied  from  reading  to  reading  were 
the  slit  width,  which  was  set  at  arbitrarily  chosen  values; 
the  dark  current  dial  on  the  DU,  which  set  the  base  line  at 
zero  on  the  chart;  the  sensitivity  setting  on  the  DU  power 
supply,  which  was  adjusted  to  keep  the  signal  on  the  chart 
paper  throughout  the  range  of  slit  widths  and  concentrations. 

To  establish  the  recorder  base  line,  the  following 
procedure  was  observed:  the  shutter  was  closed;  the  blank 

signal  with  distilled  water  only  being  aspirated  into  the 
flame  was  measured;  the  signal  (and  noise)  was  measured  by 
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allowing  the  sample  solution  to  aspirate  approximately  two 
full  minutes  (an  inch  on  the  chart  paper).  In  Figure  1 a 
typical  chart  reading  of  the  base  line,  blank,  and  signal 
is  given. 

The  translation  of  the  recorded  signals  into  values  of 
the  signal-to-noise  ratio  requires  some  explanation.  The 
signal  value  was  taken  to  be  the  distance  in  inches  between 
a line,  parallel  to  the  base  line  and  running  through  the 
center  of  the  noise  band  of  the  blank  signal,  and  a similar 
line  running  through  the  center  of  the  noise  band  of  the 
sample  signal.  The  peak-to-peak  noise  was  taken  to  be  the 
width  in  inches  of  that  solid  black  portion  of  the  noise 
band  of  the  sample  signal.  This  peak-to-peak  width  wag 
then  divided  by  two  in  order  to  consider  that  half  of  the 
noise  band  lying  above  the  signal  line  and  also  divided  by 
the  square  root  of  two  to  convert  it  to  r.m.s.  noise  (see 
Figure  1). 

The  boundaries  of  the  solid  black  portion  of  the  noise 
band  were  obtained  by  use  of  considerable  judgment  on  the 
part  of  the  operator.  In  fact,  this  step  introduced  suf- 
ficient error  into  the  value  of  the  measured  signal-to-noise 
ratio  to  partially  nullify  any  argument  that  might  be  based 
on  the  exact  agreement  between  theoretically  determined  and 
measured  signal-to-noise  versus  slit  width  curves.  The 
estimation  of  the  noise  boundaries  is  necessary  because  of 
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Peak  to  peak  noise 


Recorded  noise 
band  width 


Sample  signal 


Water  signal  (blank) 

-A 


Fig.  1. -Typical  chart  readout,  with  interpretation,  for 
measurement  of  signal  and  noise. 
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the  following  factors:  (a)  there  are,  at  some  times, 

spurious  peaks  in  the  noise  pattern  reaching  well  beyond 
the  solid  portion;  (b)  also,  at  some  times,  there  was  an 
apparently  very  low  frequency  noise  wave  pattern  which 
lifted  and  lowered  the  whole  solid  portion;  (c)  most  of  the 
time,  the  boundary  of  the  solid  band  was  not  a sharp  line 
but  irregularly  saw-toothed.  The  total  area  between  the 
signal  line  as  the  base  boundary  and  the  pen  mark  (or  noise 
line)  as  the  upper  boundary  represents  the  r.m.s.  noise 
energy.  It  was  assumed  that  the  solid  part  of  the  noise 
band  represented  the  great  majority  of  this  area  because 
the  solid  portion  was  formed  by  many  peaks  overlapping  each 
other  due  to  the  slow  chart  speed.  Consequently,  the 
spurious  peaks  above  this  line  represented  an  insignificant 
percentage  of  the  total  area  to  warrant  inclusion  of  their 
area  in  the  measurement.  The  low  frequency  wave  motion 
affecting  the  whole  noise  band  resembled,  and  was  assumed 
to  be  due  to,  the  pattern  shown  when  certain  harmonics  of 
higher  frequency  noises  are  being  passed,  by  the  pen  speed 
filter.  In  such  cases,  it  was  assumed  that  the  noise  energy 
that  fell  within  the  36  cps  band  width  range  was  still 
represented  by  the  solid  black  portion  (even  though  it  was 
undulating)  and  so  the  width  of  the  solid  black  part  of  the 
band  was  reported  as  the  peak-to-peak  noise.  The  lack  of  a 
sharp  boundary  for  the  solid  band  resulted  in  a fixed  error 
in  reading  the  width  of  this  band  in  inches. 


59 


Tiie  width,  of  the  solid  hlack  noise  band  was  measured 
visually  with  the  aid  of  a micrometer  which  read  to  the 
thousandth  of  an  inch.  The  micrometer  was  laid  on  the  chart 
paper  and  the  reader's  eye  passed  back  and  forth  from  the 
edge  of  one  micrometer  jaw  to  the  edge  of  the  other  jaw 
until  it  was  felt  that  parallax  had  been  minimized  at  both 
jaw  surfaces.  The  width  of  the  noise  band  on  each  signal 
was  measured  at  several  points  with  the  micrometer  until  a 
satisfactory  agreement  among  the  measurements  had  been 
reached. 

Attempts  were  made  to  measure  the  noise  by  other  methods 
than  that  described  above.  An  attempt  was  made  to  measure 
the  electrical  energy  of  the  noise  signal  with  a r.m.s. 
frequency  meter,  but  it  was  not  possible  to  obtain  a meter 
that  would  measure  energy  down  to  near  dc  frequencies.  An 
attempt  was  also  made  to  measure  the  total  area  between  the 
signal  boundary  and  the  upper  noise  line  boundary  by  using 
fast  chart  speed  (four  inches  per  minute)  and  a planimeter. 

It  was  found  that  the  total  area  was  too  small,  and  the 
shape  too  irregular,  to  obtain  sufficient  accuracy  to  make 
such  tedious  work  worthwhile.  The  r.m.s.  noise,  measured 
as  described  above,  was  multiplied  by  the  correction  factor 
to  compensate  for  variations  due  to  the  decade  of  the  chart 
paper  in  which  the  signal  fell  (described  in  the  section  on 
Af)  and  then  used  as  the  denominator  of  the  signal-to-noise 


ratio. 
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It  was  estimated  that  an  error  of  0.02  inches  was 
introduced  into  the  peak-to-peak  noise  measurement  of  the 
noise  hand  because  of  the  lack  of  a sharp  boundary  for  this 
portion  of  the  band.  The  percent  of  error  that  this  fixed 
reading  error  introduced  was  calculated  for  each  reading, 
and  the  corresponding  percent  error  that  it  caused  in  the 
value  of  the  signal-to-noise  ratio  was  calcuated.  These 
percentage  error  values  are  listed  along  with  the  measured 
signal-to-noise  ratio  values  in  Tables  14  through  28,  and 
the  effect  of  the  error  is  also  indicated  in  the  graphs  of 
signal-to-noise  ratio  versus  slit  width  that  are  reproduced 
in  this  work. 

C.  Comparison  of  the  Theoretical  and  Experimental  Plots 
of  "signal -t'o-Noise  Ratio  Versus  Slit  Vidt3i~~ 

The  purpose  of  this  study  was  to  determine  the  size  of 
the  optimum  slit  width  in  various  situations  as  well  as  the 
sensitivity  of  the  signal-to-noise  ratio  to  the  slit  width, 
that  is,  the  shape  of  the  maximum  of  the  signal-to-noise 
ratio  versus  slit  width  curve.  It  was  further  intended  to 
demonstrate  the  validity  of  Equation  (17),  and  therefore  the 
validity  of  further  conclusions  drawn  concerning  the  in- 
fluence of  individual  parameters,  by  showing  agreement  in 
both  shape  and  magnitude  between  the  calculated  and  measured 
plots  of  signal-to-noise  ratio  versus  slit  width.  The 


61 


pq 

Hi 


P 

I — l 

Ph 

PQ 

P 

o 

co 

&! 

O 

M 

EH 

«i 

« 

EH 

P0 

pq  pq 

O 2 
£5  S! 
O P 
O p 

pq  Bh 
pq  m 

« d 
W P 
Eh  X 
o 

« J. 

o Sd 

P Pq 

a 
m o 
o p$ 
i — i p 
EH  Xh 
' W 


o 
pq  H 
CO  « 
H EH 

O Pq 

Sa  2 

I O 

O H 
Eh  W 
I o 

33 

2 eh 

e>  co 

H 

co  -=d 

P !P 

SH 

p 

CO 

<d 


P 

*d 

P 

pq 

EH 

<1 


O 

3 


d 

0 

d* 

CM 

fft 

ift 

lft 

ft- 

lft 

cO 

a 

• a 

d 

+i 

+ 1 

+ 1 

+ 1 

+ 1 

+1 

+1 

rH 

p 

0 

r 

i 

p 

cO 

d 

to 

CD 

CO 

cO 

fft 

CM 

Lft 

Cft 

d- 

rH 

•H 

d 

•2 

d" 

d- 

d* 

Lft 

Lft 

Lft 

CO 

cO 

s 

0 

O 

fft 

i — 1 

O 

s 

0 

O d 

o 

> 

rH  0 

cO 

-P 

• 

£ 

0 

o 

1 1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

• r» 

3 

rH 

O 

o 

O 

o 

lft 

CM 

Cft 

Lft 

o 

— 1 

O 

lft 

Lft 

Lft 

lft 

Lft 

d* 

d- 

fft 

rH 

c* 

1 

rH 

ift 

d 

0 

•rH 

o 

0 

a 

p ift 

0 

a 

a 

o 

d 

• 

0) 

ft 

CM 

CM 

CM 

CM 

cO 

cO 

Ift 

s 

CM 

ift 

fH 

o 

ft- 

• 

d 

+1 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

H- 

1 — 1 

• 

o 

0 

P 

1 — 1 

0 

p 

o 

d- 

d- 

d- 

ft- 

fft 

ift 

ft- 

d 

0 

0 

♦ 2 

i — l 

rH 

CM 

CM 

lft 

Lft 

d* 

rft 

0 

P 

O 

fcO 

0 

d 

o 

o 

fH 

0 

H d 

U 

P 

0 

d 

d 

0 

P 

X 

o 

a 

0 

P 

•H 

•rH 

rH 

p 

p 

d 

IN 

o- 

CO 

CO 

d- 

Lft 

Cft 

c0 

ift 

0 

0 

o 

d* 

d" 

d* 

d" 

d" 

fft 

CM 

CM 

rH 

W 

a 

0 

rH 

0 

•H 

0 

pft- 

O 

• 

CO  CM 

ift 

0 

c0 

CO 

0 

(J 

d 

Lft 

cO 

oo 

d“ 

cO 

d* 

CM 

rH 

0 

• 

• 

• 

• 

• 

• 

p 

<■* 

a 

co 

a 

d 

+ 1 

+ 1 

+1 

+i 

+ 1 

+ 

CM 

cti 

o 

to 

i 

i 

0 

0 

p 

0 

* 

• 

• 

• 

♦ 

• 

a 

o 

0 

d" 

cO 

00 

co 

ft- 

a 

• ** 

0 

2 

P rH 

P 

• 

CO 

1 

O 

(H 

d 

h> 

rH  d 

0 

•H 

o 

0 

P 

a 

d 

P 

S Ift 

0 

d 

o 

0 

a 

•rH 

• 

p 

o 

o 

d 

o 

00* 

CD 

d- 

lft 

d" 

o 

CO 

ift 

•rH 

o 

rft 

ftJ 

CM 

CM 

CM 

rH 

rH 

rH 

0 

• 

p 

§ 

o 

p 

d 

3 

o 

0 

O 

» — i 

Lft 

P 

#* 

d 

a 

p 

CM 

o 

p 

p 

d 

•rH 

ft 

d 

0 

0 

p 

•rH 

p 

P 

0 

/^*N 

cO 

0 

N 

a x a. 

S 

a 

•H 

o d 

p 

•H 

a 

•h  a 

•rH 

p 

d 

o 

•p  -p  P 

rH 

0 

o 

p 

d co  P 

0 

pq 

i — 1 

0 

a d 

P 

0 

P 

P -H 

S'-' 

• r* 

d Os 

O 

O 

O 

O 

O 

o 

o 

o 

a a. 

d o-d 

0 

fft 

eft 

ft- 

Cft 

rH 

o 

CO 

o 

0 

O +5 

i — 1 

d* 

co 

o 

•rH 

faO 

cO 

d *h 

rH 

p 

h> 

ift 

p 

X 

Lft 

o co 

o 

o 

lft 

Calculated  optimum  slit  width 


62 


A 

rH 


no 

£ 


EH 


M 

Ph 

<4 

m 

fo 

o 

CO 

S3 

o 

M 

EH 

a 

EH 

S3 

PO 

O 

S3  0 

O no 
o 


% 

P 

ft 


no 
no 
B 

Sd  S3 
Eh  no 
t5 
W >h 
o x 

Ph 


o 

I 

S3 

no 

C5 
O 
ft 

o 

ft 

niW 

CO 

h w 

o o 

H 


CO 

o 

H 

Eh 

< 

Ph 


I 

o 

EH 

I 


Ph 


(3 

Ct> 

H 

CO 

ft 

M 

B 

p 

co 


P 

no 


S3 

H 


P 

S3 

«f 

P 

no 

EH 

<o 


o 


o 


P 

% 

• 

• 

« 

• 

(A 

• 

A 

0 

d" 

cD 

A 

A 

ft 

rH 

co 

rH 

• 0* 

fn 

+1 

— 1 

d 

+| 

+ 1 

+ 1 

+1 

+1 

+1 

+ 

i 

p 

0 

d 

■p 

0 

•H 

hO 

© 

A 

CM 

A 

A 

o 

A 

CO 

A 

a 

d 

• a 

A 

cD 

CD 

ft 

ft 

A 

ft 

co 

A 

0 

o 

a 

i — i 

o 

o 

0 

O P 

> 

i — 1 0 

• 

0 

-p 

o 

d 

0 

o 

1 — 1 

o 

d 

d" 

d" 

A 

A 

A 

A 

A 

1 — 1 

o 

H 

o 

o- 

ft 

ft 

ft 

ft 

CD 

CD 

CD 

CM 

o 

1 

rH 

H 

d 

0 

•H 

O 

0 

a 

-P  A 

0 

a 

« 

fH 

o 

• 

« 

• 

• 

• 

% 

CM 

• 

• 

0 

C\J 

A 

CM 

ft 

A 

A 

rH 

A 

s 

CD 

A 

fH 

o 

A 

• 

d 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+1 

+ 

rH 

• 

O 

m 

ft 

rH 

0 

P 

0 

d* 

CM 

CM 

CM 

CO 

i — 1 

O 

d 

0 

0 

• s 

A 

d" 

d* 

CD 

A 

CD 

ft 

A 

0 

p 

p 

o 

faO 

0 

o 

o 

fH 

0 

rH 

fH 

p 

0 

P 

d 

0 

p 

>5 

o 

a 

0 

P 

•H 

•H 

• 

-p 

p 

d 

rH 

o 

A 

CO 

ft 

A 

CM 

O 

•** 

0 

0 

o 

ft- 

ft 

ft 

CD 

CD 

A 

d“ 

d" 

CM 

M 

fH 

0 

0 

•rH 

0 

ftft 

o 

• 

0 CM 

A 

0 

d- 

Cl) 

CD 

0 

rd 

CM 

rH 

0 

CM 

rH 

CM 

A 

A 

CM 

rH 

ft 

0* 

fH 

CO 

s 

fH 

d 

+ 1 

+ 1 

+ | 

+ 1 

+ 1 

+ 1 

+ 1 

+ 

CM 

3 

o 

CO 

0 

0 

p 

0 

fH 

o 

0 

O 

CD 

CD 

rH 

cO 

A 

cD 

A 

d 

• 0* 

0 

s 

rH 

rH 

1 — 1 

CM 

CM 

rH 

rH 

rH 

-p 

rH 

ft 

« 

0 

i 

o 

U 

d 

rH 

0 

•rH 

o 

0 

ft 

a 

d 

-p 

a 

A 

0 

0 

0 

a 

•H 

• 

p 

o 

o 

d 

cD 

A 

O 

ft 

A 

A 

o 

■d" 

o 

•rH 

o 

A 

A 

A 

d* 

•d- 

CM 

CM 

rH 

CM 

0 

♦ 

ft 

a 

o 

ft 

ctf 

0 

o 

0 

o 

1 1 

A 

/^N 

ft 

0+ 

d 

s 

P 

CM 

o 

ft 

P 

P 

•H 

ft 

P 

0 

0 

p 

•H 

P 

P 

0 

d ^ 

/— N 

0 

0 

CH 

A 

a 

fH 

•rH 

o 

-p 

•H 

a 

•H  j-f 

•H 

■p 

£ 

o 

p +? 

,d 

i — 1 

© 

O 

p 

0 to 

-p 

CO 

no 

rH 

0 

fH 

p 

p 

0 

p 

-P 

B ^ 

• 0* 

d 

O 

O 

O 

o 

O 

O 

O 

o 

d d 

d 

0 

A 

A 

ft 

A 

1 — 1 

O 

O 

o 

B ^ 

0 

CD 

o 

p 

1 — 1 

d- 

co 

o 

•rH 

faO 

d 

rH 

p 

ft 

A 

o 

ft 

ri 

A 

o 

co 

o 

o 

A 

Calculated  optimum  slit  width 


TABLE  16 

CALCULATED  AND  MEASURED  SIGNAL-TO-NOISE  RATIOS  FOR  THREE  CONCENTRATIONS  OF  BARIUM 

IN  A STOICHIOMETRIC  ACETYLENE-OXYGEN  FLAMEa 
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CALCULATED  AND  MEASURED  SIGNAL-TO-NOISE  RATIOS  FOR  THREE  CONCENTRATIONS  OF  CALCIUM 
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agreement  between  curve  shapes  and  magnitude  was  not  exact 
but  certain  conclusions  can  nevertheless  be  drawn  from  a 
comparison  of  the  calculated  and  measured  values  and  curve 
shapes.  For  one  thing  the  calculated  curves  predict,  and 
the  measured  curves  bear  out  (see  Tables  21,22,23,24,25) 
that  the  signal-to-noise  ratio  for  low  concentrations  is 
dependent  on  slit  width  whereas  it  is  independent  of  slit 
width  for  high  concentrations  over  a wide  range  of  slit 
width  values.  Other  conclusions  to  be  drawn  from  a compari- 
son of  the  calculated  and  measured  curves  are  not  as 
immediately  obvious,  but  these  will  be  pointed  out  in  the 
following  paragraphs. 

If  Equation  (17)  is  used  to  study  the  effect  of  indi- 
vidual parameters,  then,  for  low  sample  concentrations, 
quite  accurate  values  for  the  two  flame  parameters,  Ic  and  If) 
must  be  known  in  order  to  predict  accurately  the  shape  of 
the  signal-to-noise  ratio  versus  slit  width  curves  (see 
Figures  5 and  6).  To  predict  the  magnitude  of  the  signal- 
to-noise  ratio,  an  accurate  value  of  Af  must  also  be  known. 
The  effect  of  I diminishes  as  the  sample  concentration 
increases,  but  the  effect  of  $ is  felt  at  all  concentra- 
tions and  is  more  significant  anyway.  In  this  study,  a 
single  value  was  assigned  to  ijJ  for  each  flame  type.  That 
value  was  based  on  average  readings  made  in  this  study. 
However,  the  complete  series  of  measurements  covered  a 
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number  of  days  and  tbe  possibility  exists  that  the  flicker 
factor  value  did  not  remain  constant  for  each  type  of  flame 
for  that  long.  If  it  is  assumed  that  ft  has  an  incorrect 
value,  then  a disagreement  would  be  expected,  in  the  shape 
and  position,  between  the  calculated  and  measured  curves  of 
signal-to-noise  ratio  plotted  against  slit  width.  The 
general  characteristics  of  the  two  sets  of  curves  would 
still  be  the  same;  they  would  all  start  at  the  origin,  rise 
to  some  maximum  height,  and  then  either  level  off  or  drop 
slightly  for  high  concentrations.  Any  deviation  from  theory 
would  fall  into  one  of  three  patterns:  (a)  for  those 

measurements  in  which  the  actual  value  of  ft  was  greater 
than  the  assumed  value  of  ft  , the  calculated  curve  would 
be  in  general  higher  than  the  measured  curve,  and  the  maxi- 
mum of  the  calculated  curve  would  be  to  the  right  (wider 
slit  widths)  of  the  maximum  of  the  measured  curve;  (b)  for 
those  measurements  in  which  the  actual  value  of  ft  is  less 
than  the  assumed  value  of  ft  , the  calculated  curve  would 
in  general  be  lower  than  the  measured  curve  and  the  maximum 
of  the  calculated  curve  would  be  to  the  left  (narrower  slit 
widths)  of  the  maximum  of  the  measured  curve;  (c)  for  those 
measurements  in  which  the  actual  and  the  assumed  values  of 
are  approximately  equal,  there  would  be  general  agreement 
between  the  two  curves.  If  only  a portion  of  the  complete 
signal-to-noise  ratio  versus  slit  width  curve  was  either 
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calculated  or  measured  (as  was  the  case  sometimes  in  this 

study)  then  the  maximum  might  fall  outside  of  the  plotted 

portion  and  the  location  of  the  maximum  would  have  to  be 

estimated  by  the  slope  of  that  portion  of  whichever  leg  of 

the  curve  appears  in  the  plot.  It  should  also  be  expected 

that  inaccuracies  in  I and  Af  values  would  alter  the 

c 

curve  shapes  to  some  extent,  but  the  curve  shapes  should 
still  be  discernable  because  the  influence  of  ^ should 
predominate. 

A complete  listing  of  the  calculated  signal-to-noise 
ratio  values  and  their  corresponding  measured  values  is 
given  in  Tables  14  through  28.  As  expected,  the  reasons 
for  disagreement  between  theory  and  experiment  are  not 
entirely  clear  in  all  cases.  Representative  plots  of 
signal-to-noise  ratio  versus  slit  width  are  given  in 
Figures  2,  3,  and  4.  In  Figure  2,  the  calculated  and 
measured  signal-to-noise  ratio  versus  slit  width  curves  for 
the  0.1  ppm  sodium  solution  in  the  stoichiometric  acetylene- 
oxygen  flame  are  compared.  This  is  considered  to  be  a case 
in  which  the  assumed  value  of  ijJ  is  smaller  than,  the  actual 
value  of  lj)  , because  the  maximum  of  the  calculated  curve  is 
to  the  right  of  the  maximum  of  the  measured  curve.  Figure 
3 compares  the  calculated  and  measured  signal-to-noise 
ratio  versus  slit  width  curves  for  the  1000  ppm  lithium 
solution  in  the  stoichiometric  acetylene-oxygen  flame.  This 
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Experimentally  measured  signal-to-noise 
ratio  value  with  magnitude  of  noise  read- 
ing error  indicated 


Fig.  2. -Comparison  of  theoretical  and  experimental  signal- 
to-noise  ratio  versus  slit  width  curves  for  a 0.1 

O 

ppm  sodium  (5890  A)  solution  in  a stoichiometric 
acetylene-oxygen  flame. 
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Fig.  3* -Comparison  of  theoretical  and  experimental  signal- 
to-noise  ratio  versus  slit  width  curves  for  a 1000 

O 

ppm  lithium  (6708  A)  solution  in  a stoichiometric 
acetylene- oxygen  flame. 


Signal-to-Noise  Ratio 


Rig.  4. -Comparison  of  theoretical  and  experimental  signal - 

to-noise  ratio  versus  slit  width  curves  for  a 10  ppm 

O 

calcium  (4227  A)  solution  in  a fuel  rich  hydrogen- 
oxygen  flame. 
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is  considered,  to  be  a case  in  which  the  assumed  value  of 
is  greater  than  the  actual  value  of  'f  , because  the  maxi- 
mum of  the  calculated  curve  is  to  the  left  of  the  maximum 
of  the  measured  curve.  Figure  4 compares  the  calculated 
and  measured  signal-to-noise  ratio  versus  slit  width  curves 
for  the  100  ppm  solution  of  calcium  in  the  fuel  rich 
hydrogen-oxygen  flame.  This  is  considered  to  be  a case  in 
which  the  assumed  and  the  actual  values  of  If  are  approxi- 
mately equal,  because  of  the  general  agreement  between  the 
two  curves.  The  majority  of  the  comparisons  seem  to  be  of 
the  type  that  indicate  an  assumed  value  of  P that  is  greater 
than  the  actual  value  of  Ip  . 

The  above  results  can  be  summarized  as  follows:  if  the 

value  of  P is  too  variable  to  be  presented  in  Equation  (17) 
as  a single  value  over  a series  of  measurements,  then 
Equation  (17)  would  indicate  a disagreement  between  calcu- 
lated and  measured  signal-to-noise  ratio  versus  slit  width 
curves.  The  agreement  between  calculated  and  measured 
signal-to-noise  ratios  indicates  that  Equation  (17)  is  valid, 

as  has  been  shown  by  other  uses  to  which  essentially  the 

25 

same  equation  has  been  put  by  Vickers.  ' Furthermore,  it 
can  be  assumed  that  the  flame  flicker  factor,  !p  , is  a 
parameter  whose  value  is  not  constant  enough  such  that  a 
single  value  can  be  assumed  for  it  throughout  a lengthy 
series  of  measurements. 


V.  INFLUENCE  ON  SIGNAL-TO-NOISE  RATIO  VERSUS  SLIT 
WIDTH  CURVES  OF  INDEPENDENTLY  VARYING  INDIVIDUAL 

PARAMETERS 

The  relative  significance  of  each  of  the  instrumental 
parameters  of  Equation  (17)  can  be  found  by  use  of  the 
computer  program  for  calculating  points  on  the  signal-to- 
noise  ratio  versus  slit  width  curve.  The  evaluation  of 
parameters  used  in  Equation  (17)  has  already  been  discussed. 
Using  these  values,  the  calculated  points  on  the  signal-to- 
noise  ratio  versus  slit  width  curve  should  produce  a curve 
of  the  same  shape  as  that  measured  experimentally,  that  is, 
a normal  curve  should  result.  Those  values  for  each  para- 
meter, which  were  considered  to  be  valid  for  sodium  in  the 
stoichiometric  acetylene-oxygen  flame,  were  arbitrarily 
considered  to  be  normal  for  this  portion  of  the  study  and 
were  used  in  Equation  (17).  To  alter  any  single  one  of 
these  values  (while  holding  the  rest  at  their  estimated 
normal  value)  should  produce  an  abnormally  shaped  curve. 

The  relative  deviation  of  the  abnormal  curve  from  the  normal 
curve  should  be  a measure  of  the  relative  significance  of 
that  parameter  to  the  value  of  the  signal-to-noise  ratio. 
Such  information  could  be  of  value  to  further  theoretical 
studies  in  flame  spectrometry. 
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The  relative  significance  of  each  parameter  might  he 
predicted  from  Equation  (17)  by  purely  algebraic  means. 
However  the  availability  of  the  computer  program,  and  the 
ease  with  which  each  parameter's  value  can  be  increased  or 
decreased,  with  a complete  set  of  curve  points  produced  for 
each  variation  of  a parameter,  makes  this  sort  of  investi- 
gation uniquely  appropriate  to  this  study.  Consequently 
each  of  the  following  parameters  (or  parameter  combinations), 
lj)  , Ic,  i^,  / and  M,  was  increased  threefold  and  decreased 
threefold  while  all  of  the  other  parameters  were  held  at 
their  estimated  normal  value.  The  parameter  Wm  was  felt  to 
be  subject  to  less  drastic  fluctuations,  and  so  it  would  be 
more  useful  to  observe  the  effect  of  a twofold  increase  and 
decrease  from  its  estimated  normal  value  of  .003  cm  (esti- 

Q 

mated  from  Herrmann  and  Alkemade  ).  The  resulting  distor- 
tions in  the  shape  of  the  signal-to-noise  ratio  versus  slit 
width  curves  for  each  parameter  are  shown  in  Figures 
7,  8,  and  9«  Several  conclusions  can  be  drawn  from  the 
shapes  of  the  curves. 

In  general,  wider  slit  widths  increase  the  signal-to- 
noise  ratio. 

At  low  concentrations  the  signal-to-noise  ratio  versus 
slit  width  curves  show  a sharper  maximum.  At  high  concen- 
trations the  signal-to-noise  ratio  reaches  a maximum  value, 
levels  off,  and  is  virtually  independent  of  slit  widths. 
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100. 


O Parameter  value  threefold  greater  than  normal 
• Parameter  value  threefold  smaller  than  normal 


Fig.  5.-Signal-to-noise  ratio  versus  slit  width  curves  for 
extreme  values  of  the  flame  flicker  factor. 


Fig.  6,-Signal-to-noise  ratio  versus  slit  width,  curves  for 
extreme  values  of  the  continuum  intensity. 
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Fig.  7« -Signal -to-noise  ratio  versus  slit  width  curves  for 
extreme  values  of  the  minimum  resolving  power  slit 
width. 
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O Parameter  value  threefold  greater  than  ^ 
normal 

• Parameter  value  threefold  smaller  than  normal 


Fig.  8,-Signal-to-noise  ratio  versus  slit  width  curves  for 
extreme  values  of  the  photosensitivity  factor  and 
the  photomultiplier  tube  gain  factor. 


Signal-to-Noise  Ratio 
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O Parameter  value  threefold  greater  than 
normal 

0 Parameter  value  threefold  smaller  than  normal 


Fig.  9.-Signal-to-noise  ratio  versus  slit  width  curves  for 
extreme  values  of  the  photoanodic  dark  current. 
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The  flame  flicker  parameter  has  the  greatest  effect  on 
the  magnitude  of  the  signal-to-noise  ratio  of  any  of  the 
parameters  studied.  It  is  the  flicker  factor  (along  with  Af) 
which  determines  the  value  of  the  high  concentration  maxi- 
mum (see  Eqn.  (40)). 

Increasing  the  / and  M factors  causes  a slight  increase 
in  signal-to-noise  ratio  that  can  he  observed  for  low  sample 
concentrations  at  slit  widths  below  the  optimum.  This 
indicates  that  increased  sensitivity  amplifies  the  signal 
somewhat  more  than  it  amplifies  the  noise,  at  least  at 
narrow  slit  widths. 

Larger  values  for  the  photoanodic  dark  current,  which 
results  in  larger  instrumental  noise,  has  the  greatest 
effect  at  narrow  slit  width  and  low  concentration. 

It  is  interesting  to  note  that  a given  increase  in  the 
photoanodic  dark  current,  i^,  will  cause  a similar  change 
in  the  curve  that  would  result  from  an  identical  decrease 
in  the  photomultiplier  tube  gain  factors, V and  M. 

In  general,  instrumental  noises  are  significant  at 
narrow  slit  widths,  and  flame  noises  are  predominant  at 
normal  or  wide  slit  widths. 

The  only  parameter  which  greatly  influences  signal-to- 
noise  ratio  as  slit  widths  increase  is  the  continuum  in- 
tensity. Opening  the  slit  wider  increases  the  continuum's 
proportion  of  the  total  intensity  entering  the  instrument. 
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This  is  because  increasing  the  slit  also  increases  the 
spectral  band  width.  Consequently  Equation  (40)  in  which 
it  is  assumed  that  the  continuum  intensity  had  become  insig- 
nificant, is  no  longer  valid,  and  the  horizontal  curve  at 
the  maximum,  which  it  predicts,  does  not  occur.  Rather 
the  equation  which  does  apply  is 

i/AiT  = I/Af1/2^  [IcRd(W+Wm)+I]  (43) 

and  the  signal-to-noise  ratio  will  decrease  with  increasing 
slit  widths. 

The  two  flame  parameters,  continuum  intensity  and  the 

flame  flicker  factor,  and  the  shot  noise  term  control  the 

value  of  the  optimum  slit  width,  that  is,  the  location  of 

the  maximum  of  the  signal-to-noise  ratio  versus  slit  width 

curve.  At  low  values  for  I its  depressing  influence  is 

c 

delayed,  and  the  maximum  occurs  at  wider  slits.  When  I 

c 

is  large,  its  depressing  influence  is  felt  at  narrower  slits 
so  the  optimum  slit  width  is  smaller.  The  flicker  factor 
controls  the  proportion  of  noise  that  accompanies  the 
continuum  intensity,  and  so,  for  a given  I , higher  values 
of  p would  have  the  same  result  that  a higher  value  of  IQ 
would  have.  Consequently  large  values  for  either  I , or  jp  , 
or  both,  tend  to  cause  the  optimum  slit  width  to  be  narrower; 
small  values  for  either  Ic,  or  ij)  , or  both,  tend  to  cause 
the  optimum  slit  width  to  be  larger. 


The  difference  in  shape  between  the  high  and  low 
concentration  signal-to-noise  ratio  versus  slit  width 
curves  is  due  to  the  relative  sizes  of  spectral  line 
intensity  and  continuum  intensity.  Continuum  intensity  is 
approximately  constant,  but  spectral  line  intensity,  of 
course,  varies  with  sample  concentration.  If  the  spectral 
line  intensity  is  high,  the  continuum  intensity  will  be 
insignificant,  and  Equation  (4-0)  will  control  the  shape  of 
the  curve  until  very  wide  slit  widths  are  reached,  at  which 
point  Equation  (4-3)  will  be  valid  and  the  curve  will  begin 
to  drop.  If  the  spectral  line  intensity  is  low,  Equation 
(4-3)  will  become  valid  at  much  narrower  slit  widths,  and 
there  will  not  be  a broad  flat  maximum  to  the  curve. 


VI.  SUMMARY  OR  CONCLUSIONS 


1.  Decreasing  flame  noises  offers  more  promise  of  im- 
proving tiie  signal-to-noise  ratio  than  decreasing  instrument 
noises.  However  narrower  frequency  response  hand  widths 
will  improve  the  ratio. 

2.  A decrease  in  the  flame  flicker  noise  results  in  the 
greatest  improvement  in  the  signal-to-noise  ratio. 

3.  The  flame  flicker  factor  does  not  have  a sufficiently 
constant  value  such  that  a single  value  for  it  could  he 
assumed,  even  in  a given  flame,  over  an  extended  period  of 
measurements.  If  values  for  this  parameter  are  required, 
they  should  he  obtained  frequently  hy  measuring  the  signal- 
to-noise  ratio  of  a high  concentration  sample. 

4.  All  other  things  being  equal,  better  signal-to-noise 
ratios  result  with  wider  slit  widths;  hut  extremely  wide 
slits  will  always  lower  the  signal-to-noise  ratio  value. 

5.  Slit  widths  are  critical  at  low  sample  concentrations 
hut  are  much  less  so  at  higher  sample  concentrations. 

6.  Instrument  noises  are  significant  at  slit  widths 
narrower  than  the  optimum  slit  whereas  flame  noises  pre- 
dominate at  slit  widths  wider  than  the  optimum  slit. 
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7.  At  slit  widths  larger  than  the  optimum  slit  width 
changes  in  sensitivity  setting,  that  is,  gain  factor  of  the 
phototube,  have  no  effect  on  the  signal-to-noise  ratio. 

8.  For  studies  in  regions  in  which  the  flame  continuum 
intensity  is  higher  narrower  optimum  slit  widths  result 
than  for  regions  of  lower  flame  continuum  intensity. 

9.  When  using  less  stable  flames  narrower  optimum  slit 
widths  result  than  with  more  stable  flames. 
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APPENDIX  1 


Fortran  Program  for  Calculating  the  Total  Partial  Pressure 
of  the  Metal  in  the  Various  Forms  Present  in  the  Flame 
Gases  as  a Function  of  Flame  Type.  Molar  Concentration  of 
the  Aspirated  Solution,  and  the  Aspiration  Rate-Efficiency 
Factor  Product. 

This  program  is  designed  to  calculate,  in  atmospheres, 

the  sum  of  the  partial  pressures  of  all  forms  of  the  metal 

element  that  are  present  in  the  flame  gases.  The  parameter, 

PT,  is  used  in  a later  program  for  calculating  the  number 

of  uncombined  metal  atoms  present  per  cnr  of  flame  gases. 

The  calculations  performed  by  the  programs  are  based  on  the 

24- 

equation  given  by  Vickers 

PT  = 3.  x 1021k(  § n298/QnT)  £0  , W 

where  k is  the  Boltzmann  constant  in  cm^  atm  °K  \ $ 

is  the  aspiration  rate  in  cm^  min  \ n29g  nT  are  tiLe 

moles  of  gases  present  before  combustion  and  the  moles  of 

gases  present  after  combustion  respectively.  (In  both 

cases,  the  moles  of  aspirated  solution  are  included  and 

the  volume  units  are  the  same.)  The  factor  Q is  the  flow 

3 o —1 

rate  of  unburned  gases  in  cur  atm  K , C is  the  atomiza- 

27 

tion  efficiency  factor  applied  to  the  aspiration  rate  and 
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C is  the  concentration  of  the  aspirated  solution  in  moles 
liter’’"'".  By  the  proper  choice  of  units  for  £>298 ’ ^ can 
made  equal  to  Q.  Then  Equation  (4-h)  reduces  to 

PT  = 0.^-lh(gfc/nT)C  (45) 

Both  of  the  terms,  <§&■  and  C,  are  independent  parameters 
while  n^  depends  upon  the  types  of  fuel  and  oxygen  and  on 
their  ratio. 

The  data  used  in  this  program  consists  of  arbitrary 
values  for  fuel  flow  rate,  oxygen  (or  air)  flow  rate, 
molar  concentration  of  aspirated  solution,  and  aspiration 
rate — efficiency  factor  product.  By  means  of  this  program 
PT  values  are  calculated  for  various  solution  concentra- 
tions being  aspirated  at  various  rates  into  flames  that 
vary  both  as  to  types  of  fuel  and  oxygen  and  as  to  ratios 
of  fuel  to  oxygen  flow  rates.  A computed  GO  TO  statement 
allows  the  program  to  select  a different  subroutine  for  each 
type  of  fuel-oxygen  combination.  A GO  TO  1 statement  makes 
it  possible  to  use  separate  decks  of  data  cards  for  differ- 
ent ranges  of  fuel  flow  rates  because  the  flow  rates  used 
for  hydrogen  are  quite  different  than  those  used  for  acety- 
lene. An  IF  statement  within  the  subroutine  allows  for  the 
selection  of  different  arithmetic  statements  for  each  of  the 
three  possible  situations  of  the  flame:  fuel  rich, 

stoichiometric,  or  oxygen  rich.  The  interpretation  of  terms 
used  in  the  program  is  as  follows: 
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FUEL ( J ) 

3 _l 

Flow  rate  of  fuel  gas  in  cur  sec 

XFUEL(J) 

3 -1 

Flow  rate  of  fuel  gas  in  cm  min 

OX(X) 

3 -1 

Flow  rate  of  oxygen  or  air  in  cur  sec 

XQX(X) 

3 . -1 

Flow  rate  of  oxygen  or  air  m cur  mm 

PHSC(L) 

3 -1 

Product  of  aspiration  rate,  in  cm  min  and 

the  efficiency  factor. 

C(M) 

Molar  concentration  of  the  aspirated  solution. 

PT(L,M) 

Total  (partial)  pressure  of  the  metal  species, 
in  all  its  forms,  present  in  the  flame  gases, 
in  atmospheres. 

JJ 

Limit  of  the  list  of  fuel  flow  rates. 

XX 

Limit  of  the  list  of  oxygen  flow  rates. 

LL 

Limit  of  the  list  of  aspiration  rate-effici- 
ency factor  product  values. 

MM 

Limit  of  the  list  of  solution  concentrations. 

JJJ 

Lower  index  limit  of  loop  involving  computed 
GO  TO. 

XXX 

Upper  index  limit  of  loop  involving  computed 
GO'  TO. 

A 

3 

Surplus  fuel  gas,  after  combustion,  in  cnr 
sec”1. 

B 

3 -1 

Surplus  oxygen,  after  combustion,  in  cm  sec 

D 

3 

Water  vapor  in  the  combustion  products  in  cur 
sec”1. 

C 

Carbon  dioxide  in  the  combustion  products  in 

3 -1 

cur  sec 

H 

Nitrogen  in  cm^  sec”1. 

XX 

3 

Sum  of  the  combustion  product  gases  xn  cm 

-1 

sec 
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The  composition  of  the  data  deck  (assuming  6 values  for  the 
flow  rates  of  both  fuel  and  oxygen,  10  values  for  PHEC,  and 
8 values  for  concentration)  are  as  follows: 


JJ 

KK 

LL 

MM 

JJJ 

KKK 

FUEL(l) 

FUEL(2) 

FUEL(3) 

FUEL(4) 

FUEL(5) 

FUEL(6) 

0X(1) 

0X(2) 

0X(3) 

0X(4) 

0X(5) 

0X(6) 

PHEC ( 1 ) 

PHEC(2) 

PHEC(3) 

PHEC(4) 

PHEC ( 5 ) 

PHEC ( 6 ) 

PHEC(7) 

PHEC(8) 

PHEC(9) 

PHEC (10) 

0(1) 

0(2) 

0(3) 

0(4) 

0(5) 

C(6) 

0(7) 

0(8) 

An  example  of  the 

way  the 

computed  values  are 

printed 

follows : 

HYDROGEN/OXYGEN  FLAME 

Fuel  Flow  Rate  3000  CC/Min 
OXYGEN  FLOW  RATE  500  CC/MIN 

LOG  CONC.  MOLES/LITER  -10  -9  -8 

ASP.  RATE  X EFF.  FAC.  CC/MIN 

0.20 

0.40 

0.60  (Pj  values) 


The  Fortran  program  itself  is  as  follows: 
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DIMENSION  FUEL ( 6 ) » OX ( 6 ) ,PHEC( 10),C(8),PT(10,8) 

I READ  INPUT  TAPE  5 , I 00 , J J , KK , LL , MM , J J J , KKK , { FUEL ( J ) » J=I 
1. JJ) , ( OX ( K ) ,K=1,KK) , (PHEC(L) , L=1,LL) , (C(M) ,M=1,NM) 
WRITE  OUTPUT  TAPE  6 f 10  I , J J , KK, LL ,MM , J J J , KKK , ( FUEL  l J ) , J 
1 = 1 1 J J ) , l OX ( K ) , K= 1 » KK ) i ( PHECI L) , L= I , LL ) , { C { M ) ,M= 1 , MM J 

100  FORMAT  (61 12/C6E12.2)) 

101  FORMAT  ( 1H2,6I12/(6E12.2) ) 

CO  bO  I = J J J » KKK 

GO  TO  ( 10,20,30,40) , I 

10  DO  18  J=1,JJ 
DO  18  K=  1 , KK 

IF  ( .5»FUEL( J)-OX(K) ) 11 ,12, 13 

11  A=0. 

B=0X(K)-.5*FUEL(J) 

D=FUEL( J) 

GO  TO  14 

12  A=0. 

8=0. 

0=FUEL( J) 

GO  TO  14 

13  A=FUEL(J)~2.*0X(K) 

8 = 0. 

D=2 . *0X ( K ) 

14  XX=A+B+D 

DO  Lb  L= 1 , LL 
00  15  M= 1, MM 

15  PT(L,M)  = {.424*PHEC(L)*C(M) ) / ( XX+ ( 22. 5»PHEC { L)  ) ) 
XFUEL=60.*FUEL(J) 

X0X=60.»CX(K) 

WRITE  OUTPUT  TAPE  6,102 

102  FORMAT  ( 1H1 , 49X , 2 1HHYDR0GEN/0XYGEN  FLAME) 

WRITE  OUTPUT  TAPE  6 , 103 ,XFUEL, XOX 

DO  16  L= 1 , LL 

WRITE  OUTPUT  TAPE  6 , 104 , PHEC ( L ) , ( PT ( L , M ) ,M=1 ,MM ) 

16  CONTINUE 
18  CONTINUE 
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GO  TC  50 

20  DO  28  J=1,JJ 
DO  28  K=1,KK 

I F ( . 5 *FUEL ( J ) — . 2*QX ( K ) ) 2 1,22, 2 3 

21  A=0. 

B=.2*QX(K)-.5*FUEL( J) 

D= . 5 *FUE L ( J ) 

H= . 8 *0X ( K ) 

GU  TO  24 

22  A=0. 

8=0  • 

D=FUEL( J) 

H=.4*0X (K) 

GO  TO  24 

23  A=FUEL( J)-.4*0X(K) 

0=0. 

D=.4*0X(K) 

H=.8*QXIK) 

24  XX=A+B+0+H 
DO  25  L=  1 » LL 
DO  25  M=1,MM 

25  PT{L,M)=(.424*PHEC(L)*C(M) ) / ( XX+ ( 22. 5*PHEC ( L ) ) ) 
XFUEL=60.»FUEL( J) 

X0X=60 . *0X ( K ) 

WRITE  OUTPUT  TAPE  6,105 
105  FURMAT  ( 1HI , 51X , 18HHYDR0GEN/AIR  FLAME) 

WRITE  OUTPUT  TAPE  6, 106 , XFUEL , XGX 
DO  26  L= 1 , LL 

WRITE  OUTPUT  TAPE  6 , 104 , PHEC ( L ) , { PT { L , M ) , M= 1 , MM ) 

26  CONTINUE 
28  CONTINUE 

GO  TO  50 

30  DO  38  J=1,JJ 
DO  38  K= 1 , KK 

IF  (FUEL! J)-.4*0X(K) >31,32,33 

31  A=0. 


B=OX (K)-2.5*FUEL( J) 

G=2.*FUEL( J) 

D=FUEL( J) 

GO  TO  34 

32  A=0. 

3 = 0. 

G=2 . *FUEL ( J } 

D=FUEL { J ) 

GO  TO  34 

33  A=FUEL( J)-.4*0X(K) 

B=0. 

G=2.*(FUEL( J)-A) 

D=FUEL ( J )-A 

34  XX=A+B+G+D 
UO  35  L=1,LL 
00  35  M=1,MM 

35  PT ( L » M ) = ( .424*PHEC(L)*C(M) ) / ( XX+ ( 22 . 5*PHEC t L ) ) ) 
XFUEL=60.*FUEL( J) 

X0X=60.*0XIK) 

WRITE  OUTPUT  TAPE  6,107 
107  FORMAT  ( IH1 , 49X , 22HACETYLENE/0XVGEN  FLAME) 

WRITE  OUTPUT  TAPE  6 , 103  , XFUEL , XUX 
DO  36  L=  1 , LL 

WRITE  OUTPUT  TAPE  6 , 104  , PHEC ( L ) , ( PT l L , M ) , M= 1 f MM ) 

36  CONTINUE 
38  CONTINUE 

GO  TO  50 

40  DO  48  J=1,JJ 
DO  48  K=1,KK 

IF  (FUEL! J)-.08*0X(K) )41,42,43 

41  A=0 . 

8=0X(K)-12.5*FUEL( J) 

G=2.*FUEL(J) 

D=FUEL( J) 

H= • 8 »0X ( K ) 

GO  TO  44 
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42  A=0. 

B=0. 

G=2.*FUEL( J) 

D=FUEL( J) 

H=.8*0X IK) 

GO  TO  44 

43  A=FUEL( J)-.08*0X(K) 

8=0. 

G=2.»(FUEL(J)-A) 

0=FU  EL ( J ) -A 
H=.8*0X(K) 

44  XX=A+B+D+G+H 
00  45  1=1, LL 
00  45  M=  1 , MM 

45  PT ( L ,M ) = { .424*PHEC( L)*C(M) ) / ( XX+ { 22 . 5*PHEC ( L ) ) ) 
XFUEL=60.«FUEL(J) 

XGX=60 . *0X { K ) 

WRITE  OUTPUT  TAPE  6,108 
108  FORMAT  ( 1H 1 , 50X , 19HACETYLENE/A I R FLAME) 

WRITE  OUTPUT  TAPE  6,106,  XFUEL , XOX 
00  46  L=  1 , LL 

WRITE  OUTPUT  TAPE  6 , 1 04 , PHEC I L ) , ( PT l L , M ) , M= 1 , MM ) 

46  CONTINUE 
48  CONTINUE 

GO  TO  50 
50  CONTINUE 
GO  TO  1 

103  FORMAT  ( 1HQ , 3X , 14HFUEL  FLOW  RA TE , 4X , 1F8 . 2 , 2X , 6HCC/M I N/ 
1 04X , 16H0XYGEN  FLOW  RATE , 2X , 1 F8 . 2 ( 2X , 6HCC/M IN/// 3X , 21HL 
10G  CONC.  MOLES/LITER, 5X, 3H-10, 9X, 10X, 2H-8, 10X,2H-7,10X 
l,2H-6,10X,2H-5,10X,2H-4,10X,2H-3//3X,28HASP.  RATE  X EF 
IF.  FAC.  CC/MIN) 

106  FORMAT  ( 1H0 , 3X , 14HF UEL  FLOW  RATE , 4X , 1F8 . 2 , 2X , 6HCC/MI N/ 
1/4X.16HAIR  FLOW  RATE . 2X , 1F8 . 2 , 2X , 6HCC/MIN///3X , 2 1HL0G 
1C0NC.  MULES/LITER, 5X,3H-10,9X, 2 H- 9, 10X,2H-8, 10X,2H-7,1 
10X,2H-6,10X,2H-5,10X,2H-4,10X,2H-3//3X,28HASP.  RATE  X 
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1EFF.  FAC.  CC/MIN) 

104  FORMAT  { IHO , 5X , IF  12 . 2 , 6X , 8E 1 2 .2 ) 
END 


APPENDIX  2 


Fortran  Programs  for  Calculating  the  Number  of  Metal  Atoms 
3 

Per  cur  of  Flame  Cases,  Integrated.  Intensity  of  the  Spectral 
Line  for  Atomic  Emission,  and  Absorption  Coefficient  as 
Functions  of  Element  and  Line,  Flame  Type,  Temperature,  and 
the  Total  Partial  Pressures  of  All  Species  of  the  Metal 
Present  in  the  Flame  Cases. 

Two  programs  were  used;  one  was  designed  for  the  alkali 

metals,  and  one  was  designed  for  the  alkaline  earths.  They 

are  very  similar  to  each  other  and  differ  only  in  those 

arithmetic  statements  dealing  with  the  expressions  for  the 

dissociation  constants  of  the  metal  chlorides.  For  alkaline 

earths,  this  expression  contains  the  partial  pressure  of  the 

chlorine  atom  as  a squared  term;  the  alkali  metal  expression 

uses  the  corresponding  term  for  the  partial  pressure  of  the 

chlorine  atom  to  the  first  power.  Both  Fortran  programs  are 

given  here  in  complete  form  hut,  since  the  programs  are 

exactly  the  same  in  all  other  respects,  only  one  discussion 

of  the  program  will  he  given. 

These  programs  are  primarily  designed  to  calculate  the 

3 

number  of  uncomhined  metal  atoms  present  per  cnr  of  the 

flame  gases.  From  that  value  the  intensity  of  emitted 

-2  -1 

light  in  watts  cm  ster  , I,  and  the  absorption  coefficient 
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of  the  flame,  k°  in  cm  ' are  calculated.  The  .former  of 
these  two  is  used  in  a later  program  of  this  v/oi'k  to  calcu- 
late the  signal,  noise,  signal-to-noise  ratio,  and  optimum 
slit  width;  the  latter  will  he  useful  if  a similar  study  to 
this  one  is  done  on  atomic  absorption  flame  spectrometry. 

The  number  of  atoms  in  the  flame  depends  upon  the 
number  introduced,  P^,,  and  on  the  dissociation  constants, 
or  ionization  constant,  of  the  species  competing  with  the 
atomic  species.  These  constants  depend  upon  the  tempera- 
ture and  on  the  partial  pressures  of  the  gases  (including 
electron  pressure)  involved  in  the  formation  of  the  com- 
peting species.  The  partial  pressures  of  concern  depend, 
in  turn,  on  the  flame  type.  Consequently  the  calculated 
values  of  N,  I,  and  k°  are  produced  by  the  programs  in 
three  dimensional  arrays  for  which  temperature,  T,  and 

flame  type  are  the  independent  variables.  The  first  value 
to  be  calculated  by  the  program  is  the  partial  pressure  of 
the  atomic  species  of  the  metal,  P^.  The  technique  used  by 
the  program  for  this  calculation  is  that  of  successive 
approximation.  Equation  (JO)  for  the  alkaline  earths  i§ 
arranged  in  the  following  way: 

2L(ICXK2K3K4)  + 4(K2K?K4Pt2)  + (K^K^Pq)  + (K^K^)  + 
(K1K2K5POH)]PM  + K1K2K3K^[(Pe')2  + 4K5Prl]1/2  - 

(46) 


K1K2K3K4(Pel  + 2PT)  . 
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Equation  (31)  for  the  alkali  metals  is  similarly  rearranged; 

2[(K1K2K5Kif)  + (K2K3K^Pt)  + (K^K^)  + (KjE^Pjj)  + 

, 2 1/2 

(KlK2K3POH)]PM  + KlK2K3K4[^Pe  > + 4K5PM] 

^KgK^CPg'  + 2Pt)  . (47) 

The  coefficient  of  the  first  term  on  the  left  of  each 
equation  is  set  equal  to  the  dummy  term,  A,  for  convenience. 
The  left  side  of  the  equation  is  set  equal  to  the  dummy 
term,  Q,  and  the  right  side  is  set  equal  to  the  dummy  term, 
R.  Thus  for  any  arbitrary  value  of  the  term  P^  if 

|l  - Q/r|  < .05  , (48) 

it  is  assumed  that  this  value  of  P^  comes  very  close  to 
satisfying  the  original  equation  (either  (46)  or  (47))  and 
is  therefore  an  acceptably  correct  value  for  P^.  When  the 
arbitrarily  chosen  value  of  P^  does  not  satisfy  the  test  of 
Equation  (48),  it  is  necessary  to  determine  whether  the 
approximation  is  too  high  or  too  low.  This  information  is 
gained  by  another  IF  statement  test.  The  R term  is  inde- 
pendent of  PM.  Consequently  if 

Q > R , (49) 

it  is  because  the  approximated  value  of  P^  is  larger  than 
it  should  be  and  vice  versa.  The  complete  procedure  used 
is  to  take  P^  » 1/10  PT  as  the  initial  approximation;  test 
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this  value  for  P^  with  Equation  (48);  if  it  fails  that  test, 

an  I?  statement  branching  device  using  the  inequality  49 

causes  the  successive  approximations  to  either  increase  or 

decrease  as  the  program  proceeds  to  an  acceptable  value  of 

p.,.  In  either  case,  the  succeeding  aoproximations  take 
n 

tenfold  jumps.  The  same  IP  statement  using  inequality  49 
determines  whether  or  not  the  tenfold  jump  has  overshot  the 
true  value  of  P.,.  When  this  is  the  case,  the  program  backs 
up  to  the  last  previous  approximated  value  0f  PM  and  proceeds 
forward  again  by  twofold  jumps.  The  same  test  for  over- 
shooting is  again  employed,  and  when  it  does  overshoot,  the 
program  again  backs  up  to  the  immediately  previous  approxi- 
mation and  proceeds  forward  by  5 per  cent  jumps.  A limit 
on  the  number  of  5 per  cent  jumps  was  placed  in  the  program 
to  prevent  some  mechanical  error  from  consuming  too  much 
computer  time.  If  that  limit  should  be  reached,  the  program 
terminates  and  prints  out  information  to  that  effect. 

The  value  of  N is  calculated  from  the  value  of  P^  by 
Equation  (33)  as  described  earlier.  The  value  of  I is 
calculated  from  the  value  of  N by  Equations  (34)  or  (35),  and 
k°  is  calculated  from  the  value  of  IT  by  Equation  (38);  both 
calculations  were  also  described  earlier.  The  program  uses 
terms  which  represent  collections  of  the  constants  involved 

in  Equations  (34,  35,  3?,  and  38).  The  calculation  of  I 

32 

involves  the  "intersection  point"  value  of  N. 


This  value 
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is  dependent;  on  the  temperature,  and  so  a series  of  hand 
calculated  alues  of  the  term  including  the  "intersection 
point"  N ure  fed  in  as  data,  one  for  each  of  the  arbitrary 
temperature  values  used.  These  pieces  of  data  are  sub- 
scripted with  the  same  symbol  as  the  temperature  data.  The 
complete  data  required  by  either  of  these  two  programs 
consists  of:  a list  of  arbitrary  values  for  temperature 

sind  PT;  the  partial  pressures,  at  each  temperature  and  for 
each  flame  type,  of  hydrogen  (H) , oxygen  (0),  hydroxide 
(OH),  and  the  electron  pressure;  the  dissociation  constants 
at  each  temperature  for  the  chloride,  monoxide,  monohydride 
and  monohydroxide;  the  ionization  constant  for  the  first 
stage  of  ionization  at  each  temperature;  values  for  the 
various  collections  of  constant  terms.  In  addition  to 
these  data,  there  are  certain  labeling  directions  for  the 
printout  required.  The  data  for  one  line  of  one  element 
constitutes  a complete  data  deck  for  these  programs.  A GO 
TO  1 st...'  ement  allows  for  the  use  of  as  many  data  decks  as 
is  desired  to  be  used  in  the  same  computer  run.  The 
interpretation  of  the  program  terms  is  as  follows: 

EHIS(J,K,I)  Log  of  integrated  intensity  of  the  spectral 

-2  -1 

line  in  watts  cm  ster  calculated  by  the 
program. 

ABS(J,K,I)  Log  of  absorption  coefficient  calculated  by 
the  program. 

3 

Log  of  the  number  of  atoms  of  metal  per  cur  of 
flame  gase_  calculated  by  the  program. 


Z(J,K,I) 
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XEMIS(J,K,I) 

Integrated,  intensity  of  the  spectral  line  in 
^ _0  —1 

watts  cm  ~ ster  calculated  by  the  program,  I. 

XABS(J,K,I) 

Absorption  coefficient  calculated  by  the  pro- 

i o -1 

gram,  k cm 

XZ(J,K,I) 

Number  of  atoms  of  metal  per  cm^  of  flame 
gases  calculated  by  the  program,  N. 

XKCL(I) 

Dissociation  constant  of  the  metal  chloride. 

XKO(I) 

Dissociation  constant  of  the  metal  monoxide. 

XKH(I) 

Dissociation  constant  of  the  metal  monohydride. 

XKOrl ' “ ) 

Dissociation  constant  of  the  metal  monohydroxide. 

XKE(I) 

Ionization  constant  of  the  metal. 

PTOT(E) 

Total  of  the  partial  pressures  of  all  the 
species  of  the  metal  present  in  the  flame 
gases,  PT 

T(I) 

Arbitrary  temperatures. 

?0(J) 

Partial  pressure  of  oxygen  (0)  in  the  flame 
gases . 

PE^  J ) 

Partial  pressure  of  hydrogen  (H)  in  the  flame 
gases . 

ECHO 

Partial  pressure  of  the  neutral  hydroxide  (OH) 
in  the  flame  gases. 

PE 

Partial  pressure  due  to  electrons  in  the  flame 
gases. 

L 

Index  limit  of  the  subscript  J which  identi- 
fies the  flame  type. 

M 

Index  limit  of  the  subscript  K which  identi- 
fies the  tcv.nl  of  the  partial  pressures  of  all 
the  species  of  the  metal  present  in  the  flame 
gases . 

N 

Index  limit  of  the  subscript  I which  identi- 
fies the  temperature. 

Ill 


PM 

Q, 

H 

A 

IKY 

SAM 

BAM 

DAM 

SEEK 

XKL 

*c p\rn 

ZSA(I) 

ELEMT 

XLIME 

XFTP(J) 


Lummy  term  for  the  partial  pressure  of  the 
metal  atom  in  the  flame  gases. 

Dummy  term  for  the  value  of  the  left  side  of 
Equations  (46)  and  (47). 

Dummy  term  for  the  value  of  the  right  side  of 
Equations  (4-6)  and  (47). 

Dummy  term  for  the  coefficient  of  P^  in  the 
left  side  of  Equations  (46)  and  (47). 

Dummy  term  to  control  the  number  of  cycles  of 
approximation  used  by  the  program. 

Dummy  term  used  in  obtaining  the  antilog  of 
EMIS(J,K,K). 

Dummy  term  used  in  obtaining  the  antilog  of 
Z(J,K,I). 

Dummy  term  used  in  obtaining  the  antilog  of 
ABS( J ,K,I) . 


' 2-fcT?X02su  At£M21/2 

Log  log  L -s-p-N ) 

7.15x10  / ^^871  B^ 


;Og 


hi) 

10g[ Tj 

\ ( 


a*  A 
o bu  t 


10'  4m  B 


(T) 


Excitation  energy  in  volts. 


16m  a 7.15xlO“7  P1/2B^s 

Log  logC 2 

m21/2  t5XT  su  At  L Ao2 

Symbol  of  element  used  to  label  pages  of  output. 

Specural  line  in  angstroms  used  to  label  pages 
of  output. 

Flame  type,  used  to  label  pages  of  output. 
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RATIO (J ) 

Ratio  of  fuel  gas 
pages  of  output. 

to  oxygen, 

used  to  label 

TLMP(I) 

Temperature  values,  used  to 
in  the  output . 

label  the  tables 

The  comp 

osition  of  a complete 

data  deck  (assuming  8 tempera' 

tures,  8 

values,  and  6 flame  types)  is 

as  follows: 

L M 

N 

XKCL(l) 

XKCL ( 2 ) 

XKCL ( 3 ) 

XKCL (4) 

XKCL (5) 

XKCL (6) 

XKCL ( 7 ) 

XKCL (8) 

XKO(l) 

XK0(2) 

XK0(3) 

XK0(4) 

XKO(5) 

XK0(6) 

XKO ( 7 ) 

XK0(8) 

XKH(l) 

XKH(2) 

XKH(3) 

XKH(A-) 

me  5) 

XKH(6) 

XKH(7) 

XKH(8) 

XKOH(l) 

XK0H(2) 

XK0H(3) 

XKOH(h) 

XK0H(5) 

XK0H(6) 

XKOH( 7) 

XK0H(8) 

XKS(l) 

XXE(2) 

XKS(3) 

XKE(4) 

XKS(5) 

XKS(6) 

XKE(7) 

XKE(8) 

PTOT(l) 

PT0T(2) 

PT0T(3) 

PT0T(4) 

PT0T(5) 

PT0T(6) 

PT0T(7) 

PT0T(8) 

T(l) 

T(2) 

T(3) 

T(^) 

TC5) 

T(6) 

T(7) 

T(8) 

PO(l) 

P0(2) 

P0(3) 

P0(4) 

PO(5) 

P0(6.) 

PH(1) 

PH(2) 

PH(3) 

PH(4) 

PH(5) 

PH(6) 

I?  Oix 

P0H(2) 

POH( 3) 

P0H(4) 

PCZ(5) 

P0H(6) 

PE 

SPEK 

ZSA(l) 

ZSA(2) 

ZSA(3) 

ZSA(4) 

ZSA(5) 

ZSA(6) 

ZSA(7) 

ZSA(8) 

XXL 

EXC 

EL3HT  XLINE 

TEM?(1) 

T3MP(2)  TEMP(3) 

TEMP(^)  TEMP(5) 

TEMP(6)  TEMP(7)  TEMP(8)  XFTP(l)  XFTP(2)  XFTP(3) 

XFTP(4)  XFTP(5)  XFTP(6)  RATIO(l)  RATI0(2)  RATIO(3) 
RATIO(4)  RATIO(5)  RATI0(6) 

An  examp?,  s of  the  way  the  computed  values  are  printed 
out  is  as  follows: 

NA  5390  EX/AR  = 2/5 

ABSORPTION  COEFFICIENT 

LOG-  OF  P TOTAL  -10  -9  -8  

TEMPERATURE 

1800 

2000 

2200  (values  of  absorption  coefficient) 


The  sa  0 pattern  is  used  for  printing  out  the  values  of 
numbers  of  atoms  and  intensity.  There  is  one  such  table 
for  a single  flame  type  on  each  page. 

The  two  Fortran  programs,  for  the  alkaline  earths  and 


the  alkali  metals,  are  as  follows: 


PROGRAM  FOR  ALKAL I NE  EARTH  METALS 


Li  MENS  I ON  XKCLt 10) ,XKO( 10) , XKH( 10) , XKGHt 10),XKE(10),PT 
LOT ( 10) ,T( 10)  ,PQ{ 15) , PHI  15) , POHt 15 ) ,Z ( 15 , 10, 10) »cMIS( 15 

2.10.10)  , ZSA l 10) , ADS i 15, 10, 10 ) , XEMI SI  15 , 10, 10 ) , XABS( 1 5, 

310. 10)  ,XZ(15 ,10, 10) , TEMPI  16) ,XFTP( 15) ,RATIO( 15) 

1 READ  INPUT  TAPE  5,6, L»M,N» ( XKCL ( I ) , 1=1, N) , (XKO( I ) , 1=1, 
IN) , XKH ( I ) , 1=1,N) • ( XKCH [ I ) , I = 1 , N ) , { XKE ( I ) , I = 1 , N ) , P TOT { K 
2),K=1,M),(T(I),I=1,N),{P0(J),J=1,L),(PHIJ) , J=1,L) , (POH 
31 J) , J=1 , L) ,PE,SPEK, (ZSAt  P I = 1 , N ) , XKL , EXC 
READ  INPUT  TAPE  5,80,£Ll  .(LINE,  (TEMPI  I ) , 1 = 1, N) 

80  FORMAT  ( 2X , 2A3 , 8F6 . 0 ) 

READ  INPUT  TAPE  5 , 1 CO , ( XF TP ( J ) , J = 1 , L ) , { RAT  1 0 ( J ) , J = 1 , L ) 

100  FORMAT  ( 2 X , 1 2 A 5 ) 

WRITE  OUTPUT  TAPE  6,7,L»M»N, (XKCLl I ) , I = 1 , N ) , (XKQ{ I) , 1=1, 
IN) , XKH ( 1 ) , I=1,N) , (XKOHI I ) , 1= 1 , N ) , { XKE ( I ) , I =1 , N ) , PTOT { K 
2),K=1,M),{T(1M  = 1,N),(PC(J),J  = 1,L),(PH(J),J=1,L)  , (POH 
- i J) , J=1,L) ,PE,SPEK, (ZSACI) , I =1 , N ) , XKL , EXC 
WRIiE  OUTPUT  TAPE  6 , 85 , ELEM T , XL  I NE , l TEMP ( I ) , I = 1 , N ) 

85  FORMAT  ( 1H0 , 2A5 , 8F6 .0 ) 

WRITE  OUTPUT  TAPE  6 , 101 , (XF TP( J ) , J = 1 ,L ) , (RATIO ( J ) , J=1  , 
1L) 

101  FORMAT  (1H0.L2A5) 

DO  50  J= 1 , L 

DO  50  K=1,M 
DO  30  1=1, N 
IKY=  1 

PM=0.10*PTGT(K) 

A=2 . 0* ( XKCL ( I ) *XKG(  I )*XKH< I ) *XKOH( I ) +4 . G*XKG ( I ) * XKH ( I ) 
1*XK0H( I ) *PTOT ( K ) +XKCL ( I )*XKH( I )*XKOH( I ) *P0 ( J ) +XKCL ( I )* 
2XK0  ( I )*XK0H(  I ) *PH  ( J ) +XKCL  ( I ) »XKO(  I ) «XKH{ I ) *PQH  { J ) ) 
Ci=A*PM+XKCL(I)*XKO(  I ) *XKH { I ) *XKGH{  I)*(  ( P£*PE+4 . GO*XKE ( 
il ) *PM)**0.5) 

R=XKCL( I )*XKC  . I ) *XKH  ( I ) »XKOH(  I ) * ( P.E  +2 . 00*PTOT  ( K ) ) 

IF  (A3SF(1.00-Q/R)-G.05)23,23,3 
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3 IF  (C-R)5,5,26 
5 PM=10. 00*PK 

Q=A*PM  + XKCL ( I ) *XKG ( I )*XKHt I ) *XKGH( I )*( { PE*P£+4.CG*XK£{ 
1 1 )*PM) **0.5) 

IF  (ABSF11.00-Q/R)-0.05)23,23,9 
.9  IF  l Q-R ) 5, 1 1 , 1 1 
11  PM=0.10*PM 
13  PM=2.0G*PK 

Q=A*PM+XKCL( I )*XK0( I)*XKH( I)*XKGHCI)*( ( P£*P£+4 . 00*XKE ( 
II ) *PM) **0.5) 

IF  l ABSFtl.OO-Q/R 1-0.05)23, 23, 15 
15  IF  (G-R)  13,17,17 
17  PM=0.50*PM 
19  PM=i.05*PM 

G=.’  PM+XKCL  ( I D*XKU  ( I ) *XKH  ( I ) *XKOH(  I )*(  { PE*PE+4  . OC*XKE  ( 
1 1 ) *PM)**0. 5) 

IF  lABSFll.GC-Q/RJ-O. 05)23*23, 21 

21  I KY- I KY+ 1 

IF  ( IKY-20) 19,22,22 

22  'aRITE  OUTPUT  TAPE  6 , 10,  PM,  Q,  R,  I ,K,  J 
GO  TC  35 

26  PM=0 . 10*PM 

Q=A*PM+XKCL l ID*XKO( I ) *XKH ( I ) *XKGH( I )*( ( PE* PE+4 . 0 0 *XK E { 
1 1 ) *PM ) **0 . 5 ) 

IF  ( AOSF l 1 .OO-Q/R )-Q .05 ) 23 , 23, 27 

27  IF  (R-Q)26,28,23 

28  PM=10.00*PM 

29  PM=0 . 50  * PM 

Q-A*  PM  + XKCL ( ID*XKO( I ) *XKH(  I ) *XKGH( I )*( ( PE*PC+4 . 00*XKE ( 
II  )*PM) **0.5) 

IF  ( ABSFC 1. OO-Q/R) -0.05 ) 23, 23,30 

30  IF (R-Q) 29, 3 It  31 

31  PM=2.00*PM 

32  PM=0 .95*PM 

Q= A* PM+XKCL ( IO*XKG( I ) *XKH ( I ) *XKUH { I ) *( ( PE* PE +4 . GC*XKE ( 
1 1 ) *PM)  **0.5  ) 
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IF  ( ACSFI 1.0Q-Q/R)-0.05 ) 23, 23, 33 

33  I KY= I K Y+ 1 

IF  ! IKY-20)32,34,34 

34  kRITE  OUTPUT  TAPE  6 , 12 , PM , Q , R , I , K , J 
GO  TO  33 


23  GO  TO  35 

35  Z! J,K, I )=0.435*L0GF (PM/ ( ( 1.38E-22)*T( I ) ) ) 

A3S ( J|Kt I )=5PEK  + Z! J,K, I )-0 . 2 17*L0GF I T ( I ) ) 

IF!Z!  J,K, I )-ZSA( I ) ) 40, 40,45 
40  ErtIS(J,K»I)=XKL+Z!J,K,I  )- ( 5 . 05E3  ) *EXC/T l I) 

GO  TO  50 

45  cMIS!J,K,I)=XKL+0.5*ZSA!I ) +0 . 5*Z ! J , K , I ) - ( 5 . 0 5E3 ) *EXC/T 
i II) 

50  CONTINUE 


DO  105  J = 1 , L 
DO  105  1=1, N 
DO  105  K= 1 , M 
SAM=EMI S ( J » K , I ) 


XEMI S ! J » K, 1 ) =10. **SAM 
BAM=Z ( J , K, I i 


XZ(J#K,I )=10.**BAM 
DAM= AES { J , K , I) 


XABS(J,KfI)=10.**DAM 


105  CONTINUE 


00  113 

J = 1,L 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

RITE 

GUT PUT 

TAPE 

00  107 

I — 1 , N 

WRITE 

OUTPUT 

TAPE 

CONI INUE 

WRI  IE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

00  109 

1 = 1, N 

WRITE 

OUTPUT 

TAPE 

6, 114,ELEMT,XLINE,XFTP{ J), RATIO! J) 
6,115 
.6,118 


6 . 1 19 ,  TEMP ( I ) , ( XZ { J , K , I ) , K= 1 , M ) 

6,114, £LEMT»XLINE,XFTP(J) , RATIO! J ) 
6,116 
6,118 

6. 119,  TEMP! I ) , (XEMIS! J,K,I ) ,K=1,M) 
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i c 


in 

113 


CONTINUE 
WRITE  OUTPUT 
WRITE  OUTPUT 
WRITE  CUTPUT 
CO  111  I = 1 » N 
WRITE  GUTPUT 
CONTINUE 
CONTINUE 


TAPE  6, 114,ELEMT,XLINE,XFTP(  J)  .RATIO!  J) 
TAPE  6,117 
TAPE  6,113 

TAPE  6, 1 19, TEMPI  I ) , (XAES( J , K, I ) »K=i,M) 


6 

7 

10 

12 


GO  TO  1 

FORMAT  l 31 10/ ( 4E15. 3 ) ) 

FORMAT  (1H0,3I10/(4E15.3) ) 

FORMAT  (1X,11HPM  IS  SMALL , 3E 12 . 3 , 3 1 6 ) 
FORMAT  ( IX , l LHPM  IS  L ARGE , 3 E 12 . 3 , 3 I 6 ) 


114  FORMAT  ( 1H  1 , 4X , 4A5 ) 

115  FORMAT  { 1H0,40X,35HNUMBER  OF  ATOMS  PER  CC  OF  FLAME  GAS 


1 ) 

116  FORMAT 

117  FORMAT 
113  FORMAT 


( 1H0,42,30HINTENSITY  IN  WATTS/SQ  CM,STER.) 

( 1H0,46X,22HABS0RPTI0N  COEFFICIENT) 

{ 1H0.5X, 14HLGG  OF  P TOTAL , 8X , 3H- 10 , 10X , 2H-9 , 1 OX 
1 ,2H-8, iQX,2H-7, 1GX, 2H-6, 10X, 2H-5, 10X , 2H-410X, 2H-3//6X , 
2 1 IHTEMPERATURE ) 

119  FORMAT  C1H0,8X,1F5.0,10X,8E12.2) 


END 


PROGRAM  FOR  ALKALI  METALS 


DIMENSION  XKCLf 10) ,XKO( 1C) , XKH{ 10) ,XKOH( 10) , XKE( 10) ,PT 
lUT(lO) ,T( 10) ,P0(15) ,PH( 15) ,POH( 15),Z(15,10,10),EMIS(15 
2,10,  10)  ,ZSA( 10) ,ABS( 15, 10, 10) .XEMISl 15, 10, 10),XABS(  15, 

310, 10) ,XZ( 15, 10, 10) , TEMPI  16) ,XFTP( 15)  , RATIO!  15) 

1 READ  INPUT  TAPE  5,6,L,M,N, ( XKCL II),I=1»N),( XKOI I ) , 1=1 , 
IN) , XKH ( I ) ,I»1,N) , (XKOHt I ) , I=1,N), (XKEII ),I=1,N) ,PTOT(K 
2) ,K=1,M) , { T( 1 ) , I=1,N) , ( POl J ) , J=1,L ) , (PH( J) , J=1,L) , l POH 


118 


3( J) , J=1,L) ,P£,SPEK, { ZSA ( I ) , I = 1 , N ) , XKL , EXC 
READ  INPUT  TAPE  5 , 80 , EL EMT , XL INE , ( TEMP ( I ) , I =1 , N I 
SO  FORMAT  ( 2X , 2 A5 , 3F6 . 0 ) 

READ  INPUT  TAPE  5 , 100, ( XFTP ( J ), J=1 , L ),( RATIO ( J ), J=1 , L) 
LOG  FORMAT  (2X.12A5) 

WRITE  OUTPUT  TAPE  6 ,7,L ,M,N, IXKCLII) , 1=1 ,N ) , l XKO( I ) , 1=1 • 
IN) ,XKH( I ) , I = 1,N) , { X KGH ( I ) , I = 1 , N ) , ( XKE ( I ) , 1 = 1 , N ) ,PTOT(K 
2) ,K=1,M) , (Tl I ) , I=1,N) , (PG( J) , J=1,L) , (PH( J) , J=1,L) ,(POH 
31JI , J=1,L) i PE  » SPEK » (ZSA (I ) , 1=1 ,N ) , XKL, EXC 
WRITE  OUTPUT  TAPE  6 , 85 , ELEMT , XL INE, ( TEMP ( I ) , I = 1 ,N ) 

85  FORMAT  { 1HG , 2A5 , 8F6 . 0 } 

WRITE  OUTPUT  TAPE  6 , 10 1 , ( XFTP ( J ) , J=1 , L ) , ( RAT 10 ( J ) , J=1 , 
1L) 

101  FORMAT  ( 1H0, 12A5 ) 

DO  50  J = 1 » L 
DO  50  K=1,M 
DO  50  I = 1,N 
IKY=  1 

PM=0. 10*PTGT (K) 

A=2 . 0* ( XKCL ( i ) *XKO { I ) *XKH( I ) *XKGH( I ) + XKQ ( I ) *XKHl I ) 
i*XKUH( I ) *PTOT { K ) +XKCL ( I ) *XKH ( I > *XKQH { I ) *PG ( J ) +XKCL ( I )* 
2XK0 ( I ) *XKOH(  II *PH( J J +XKCLI  I )*XKC( I )*XKH( I ) *PGH( J) ) 
G=A*PM+XKCL( I J *XKO( I ) *XKH( I) *XKCH( I ) *( ( PE*PE+4. 00*XKE ( 

II ) *PM) **0.5 ) 

R=XKCL( I )*XKGi I ) * X K H ( I)*XKGH{ I ) * ( P E+2 . 00*P TO T ( K ) ) 

IF  (ABSFl 1 .00-Q/R5-G.05  )23,23,3 
3 IF  ( Q-R) 5,5,26 
5 PN=1C.OO*PM 

G=A*PM+XKCL ( I ) *XKO( I ) *XKH ( I ) *XKCH ( I )*{ ( PE*PE+4 . OC*XKE ( 

1 1 1 *PM ) **0. 5 ) 

IF  ( A3SF ( 1.00-Q/R) -0.05)23, 23, 9 
9 IF  ( Q-R 15,11,11 
11  PM=G . 10  * PM 
13  PM=2.00*PM 

Q=A*PM+XKCLl i ) *XKO ( I ) *XKH( I )*XKOH( I )*( ( PE*PE+4. 00*XKE ( 

1 1 ) *PM) **0.5 ) 
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1 F (ABSFll.OO-Q/RJ-O. 05)23, 23, 15 
15  IF  (Q-R)  13, 17, 17 
17  PM=G . 50* PM 
19  PM=1.05*PM 

Q=A*PM+XKCL( 1 D#XKQ ( I ) *XKH( I )»XKCH( I)M ( PE*PE+4. 00*XKE ( 

il j *PH) #*0,5 ) 

IF  i ABSF( 1.00-G/R)-0.05 )23,23,21 

21  IKY= IKY+ 1 

IF  (IKY-20) 19,22,22 

22  WRITE  OUTPUT  TAPE  6 , 10 , PP , Q , R , I , K , J 
GO  TO  35 

26  PM=0 . 10*PM 

Q=A°  PM+XKCL ( 1 D*XKQ ( I)*XKH( I )*XKOH( I)  *{ ( PE*PE+4. 00»XKE { 
1I)*PK)**0.5) 

IF  (ABSF{ 1.00-Q/R)-0.05 )23, 23,27 

27  IF  (R-Q) 26, 28,23 
2 B PM=10.00*PM 

29  PM=0.50«PM 

Q=A*PM+XKCL( ID*XKO( I ) *XKH { I ) *XKGH { I ) * ( { PE*PE+4.C0*XKE { 
1 1 ) *PM)**0.5 ) 

IF  ( ABSF ( 1 . OO-Q/R )-0. 05)23, 23, 30 

30  IF(R-Q)29,31,31 

31  PM=2.00*PM 

32  PM=0.95*PM 

Q=A*PM+XKCL( I D*XKO { I ) *XKH ( I )*XKGH( I)*( ( PE*PE+4 . GO*XKE ( 
II ) *PM) **C.5 ) 

IF  (ABSF(1.00-Q/R)-0.05)23,23,33 

33  IKY=IKY+1 

IF  { IKY-20)32,34,34 

34  WRITE  OUTPUT  TAPE  6 , 12 , PP , Q , R , I , K, J 
GO  TO  35 

23  GO  TO  35 

35  Z( J,K,I)=0.435*L0GF(PM/ ( { 1 .38E-22) *T ( I ) ) ) 

ASS(  J,  K,  I )=SPEK+Z{ J,K,  I )-0 . 2 17#  LOGF  ( T { I ) ) 

IF  { Z ( J , K , I ) — Z S A ( I) )40, 40,45 

40  LMIS  (J»K,  I )=XK.L+Z(J  ,K,  I )- ( 5 . 05E  3 ) *EXC  / T { I ) 
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GO  rc  5C 

45  EMISI J,K, I )=XKL+0.5»ZSA( I)+0.5*Z( J,K, I ) - ( 5 . 05E 3 ) *EXC / T 
id) 

50  CONTINUE 

DU  105  J=1,L 
DU  105  1 = 1, N 
DO  i 05  K=1 , M 
SAM=EMIS( J,K, I ) 

XEMISC J,K,I)=10.**SAM 
BAM=Z( J,K« I ) 

XZ( J,K, I ) =10 . * * BAM 
OAM= ABS { J, K , I ) 

XABS(J,K»1)=10***DAM 
105  CONTINUE 

DO  113  J=  1 » L 

WRITE  OUTPUT  TAPE  6 , 1 14 , ELEMT , XL INE , XFTP { J ) , RAT 10 { J ) 
WRITE  GUTPUT  TAPE  6,115 
WRITE  OUiPUT  TAPE  6,118 
DO  107  I = 1 , N 

WRITE  OUTPUT  TAPE  6 , 1 19 , TEMP { I ) , ( XZ ( J , K , I ) , K=1 , M ) 

107  CONTINUE 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRI TE  OUTPUT  TAPE  6,118 
DO  109  I = 1 , N 

WRITE  GUTPUT  TAPE  6 , 119 , TEMP ( I ) , ( XEM IS ( J»K , I } , K=1 »M) 
109  CONTINUE 

WRITE  OUTPUT  TAPE  6 , 1 14 , ELEMT , XL INE , XFTP ( J ) , RAT  1 0 ( J ) 
WRITE  OUTPUT  TAPE  6,117 
WRITE  OUTPUT  TAPE  6,118 
DO  111  1=1, N 

WRITE  OUTPUT  TAPE  6 , 1 1 9 , TEMP ( I ) , ( XABS ( J , K , I ) , K=i , K ) 

111  CONTINUE 
113  CONTINUE 
GO  TO  1 

6 FORMAT  ( 3 1 10/ ( 4E15. 3 ) ) ' 


6, 114, ELEMT, XL INE , XF  TP ( J ) , RATI  0 ( J ) 
6,116 


121 


7 FORMAT 
10  FORMAT 
12  FORMAT 

114  FORMAT 

115  FORMAT 


( 1H0.3I 10/ (4E15.3) ) 

(lX.liHPM  IS  SMALL, 3EI2. 3,316) 

( lXfllHPM  IS  LARGE,3E12.3,3I6) 

( 1H1 , 4X, 4A5 ) 

11H0.40X.35HNUMSER  OF  ATOMS  HER  CC  OF  FLAME  GAS 


i) 

116  FORMAT  l 1H0, 42, 30H INTENSITY  IN  'XATTS/SQ  CM , STER  . ) 

117  FORMAT  { 1H0,46X,22HABS0RPTI0N  COEFFICIENT) 

113  FORMAT  ( 1H0,*>X,  14HLGG  OF  P TOTAL  ,8X  , 3H-10,  10X.  2H-9,  10X 
l,2H-8,10X,2H-7, 10X.2B-6, 10X , 2H-5 , 10X, 2H-4 10X , 2H-3//6X, 
211HTEMPERATURE) 

119  FORMAT  ( IH0,8X, 1F5.0, 10X,8E12.2) 

END 


APPENDIX  3 


Fortran  Program  for  Calculating  the  Signal,  Noise,  and 
Signal-to-Noise  Ratio  in  Atomic  Emission  Flame  Spectrometry 
as  Functions  of  Element  and  Line,  Flame  Type,  Temperature, 
Slit  Width,  and  the  Total  Partial  Pressure  of  All  Species 
of  the  Metal  Present  in  the  Flame  Gases;  and  for  Calculating 
the  Optimum  Slit  Width  as  a Function  of  Element  and  Line, 
Flame  Type,  Temperature,  and  the  Total  Partial  Pressure  of 
All  Species  of  the  Metal  Present  in  the  Flame  Gases. 

This  program  uses  the  integrated  intensity  values,  I, 
produced  by  the  previous  program,  and  it  produces,  for  a 
given  element  and  spectral  line,  two  different  types  of 
calculated  values.  The  first  type  is  the  optimum  slit 
width  which  is  a function  of  flame  type,  temperature,  and 
P^,  The  other  type  includes  three  calculated  values,  signal, 
noise,  and  signal-to-noise  ratio,  and  is  a function  of  four 
different  variables:  flame  type,  temperature,  P^,  and  slit 

width. 

The  calculation  of  the  optimum  slit  width  uses 
Equation  (18)  rearranged  to  include  substitutions  according 
to  Equations  (12),  (13),  and  (14)  is  as  follows:. 
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2*fKo 2f{  + + (W»u5^  ' 

Kdaf(  W<Wm>*  I]  + 2id  ] • <5°> 

The  program  employs  a series  of  successive  approximations 
for  the  value  of  the  optimum  slit  width,  WQ.  A testing 
device  similar  to  that  used  in  the  previous  program  is 
employed  to  determine  when  an  acceptable  value  of  WQ  is 
reached.  The  left  side  of  Equation  (50)  is  set  equal  to 
the  dummy  term,  A,  and  the  right  side  is  set  equal  to  the 
au  my  term,  B.  Consequently  when  an  approximated  value  of 
V.  .satisfies  the  equation 

jl  - A/B  | < .05  , (51) 

then  that  value  of  VQ  is  considered  an  acceptable  value  for 
V, . The  initial  approximation  for  VQ  is  deliberately  set 
low,  10  microns,  and  in  the  series  the  approximation  is 
increased  by  small  increments.  The  optimum  slit  width 
should  correspond  to  the  maximum  in  the  curve  of  signal-to- 
noise  ratio  plotted  against  slit  width.  The  maximum  of  such 
a curve  is  sometimes  broad  and  flat.  This  approximation 
method  for  calculating  the  optimum  slit  width  will  only  give 
the  maximum  of  that  curve  and  will  not  give  any  indication 
of  how  broad  or  how  narrow  the  curve  is.  The  optimum  slit 
width,  W , appears  in  term  A raised  to  the  third  and  fourth 
powers  and  so  a small  change  in  the  WQ  value  used  in  the 
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series  of  approximations  is  magnified  in  its  effect  on  the 
rate  of  approach  toward  a VQ  value  that  will  satisfy 
Equation  (50).  This  rapid  approach  brings  on  the  possi- 
bility of  overshooting  the  acceptable  values  of  VQ  unless 
the  increments  of  increase  are  extremely  small.  However, 
extremely  small  increments  will  often  consume  a large  amount 
of  computer  time.  Therefore  a moderately  small  increment  is 
used,  5 microns,  with  a testing  device  to  tell  when  the 
approximated  value  of  WQ  does  overshoot  the  acceptable 
range.  The  testing  device  is  similar  to  that  used  for  the 
same  purpose  in  the  previous  program.  A low  approximation 
for  will  always  cause  A to  be  smaller  than  B because  WQ 
is  less  than  one  cm  and  is  raised  to  higher  powers  in  A than 
it  is  in  B.  Because  this  approximation  series  starts  low, 

A will  be  less  than  B at  the  beginning;  if  some  approxi- 
mated value  of  VQ  is  subsequently  used  which  makes  A larger 
than  B,  then  the  approximation  series  has  overshot  the 
acceptable  range  for  V . The  testing  device,  therefore,  is 
simply  to  test  whether  or  not  A is  larger  than  B for  every 
step  of  the  approximation.  Whenever  the  test  shows  A to  be 
larger  than  3,  the  approximation  series  terminates  and  a 
separate  statement  is  printed  out.  In  other  words,  for 
this  particular  set  of  values  of  the  variables,  the  WQ  value 
printed  in  the  list  of  WQ  values  exceeds  the  acceptable 
value  by  less  than  5 microns  (the  size  of  a single  increase 
increment) . 
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The  calculation  of  the  signal,  noise,  and  signal-to- 
noise  ratio  involves  no  approximation  and  is  simply  a matter 
of  evaluating  Equations  (15),  (16),  and  (17)  when  a particu- 
lar set  of  values  for  all  of  the  terms  is  substituted  into 
it. 

The  arrangement  of  the  data  must  provide  for  including 
four  different  variables  when  a maximum  of  three  subscripts 
per  term  is  allowed.  This  is  done  by  using  separate  data 
decks  for  each  value  of  the  fourth  variable,  the  flame 
type.  The  data  deck  treatment  can  be  somewhat  simplified 
because  some  of  the  values  to  be  used  in  the  data  deck  are 
instrument  parameters  or  properties  of  the  element  and  line 
and  are  independent  of  flame  type.  It  was  therefore  possible 
to  use  two  different  types  of  data  decks  in  this  program; 
one  which  contains  only  values  that  remain  constant  for  a 
given  element  and  line,  and  one  which  is  controlled  by  a 
loop  and  which  contains  values  that  vary  with  each  flame 
type  as  well  as  with  element  and  line.  A GO  TO  1 statement 
makes  it  possible  to  use  several  complete  sets  of  decks, 
one  for  each  element  and  line,  in  the  same  computer  run. 

The  data  used  includes  the  same  arbitrary  values  of  tempera- 
ture, P^,  and  flame  type  used  earlier  and  also  an  arbitrary 
set  of  slit  widths.  It  also  includes  the  values  of  terms 
representing  collections  of  constants  relating  to  instru- 
ment, element,  line,  and  flame  properties,  as  well  as  labels 

and  directions  for  printing  out  the  calculated  values. 
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The  scheme  used  in  printing  out  the  calculated  values, 
so  as  to  show  the  effect  of  all  variables  conveniently, 
places  each  different  flame  type  on  a separate  page.  The 
optimum  slit  width  values  then  have  only  two  variables  re- 
maining which  are  placed  in  a two  dimension  table  on  each 
page.  The  signal,  noise,  and  signal-to-noise  ratio  have 
three  remaining  variables  after  the  flame  types  are  placed 
on  separate  pages.  These  are  arranged  in  tables  with  two 
different  levels  of  horizontal  indexing  and  appear  as  one 
complete  table  per  page. 

The  interpretation  of  the  terms  used  in  the  program 


is  as  follows: 

XI(J) 

-2  -1 

Intensity  values  in  watts  cm  ster 

¥0(J) 

Optimum  slit  width  in  cm  calculated  by  the 
program. 

XW(K) 

Arbitrary  slit  width  values,  in  cm,  to  be 
used  in  calculating  the  signal,  noise,  and 
signal-to-noise  ratio. 

SIG( J,K) 

Signal  values,  in  amperes,  calculated  by  the 
program. 

RAC(J,K) 

Noise  values,  in  amperes,  calculated  by  the 
program. 

SNR(J,K) 

Signal-to-noise  ratio  values  calculated  by 
the  program. 

SW(K) 

Slit  width  values,  in  microns,  used  in  labeling 
the  tables  of  signal,  noise,  and  signal-to- 
noise  ratio. 

D(J) 

Log  of  Pm  values  used  in  labeling  the  tables 
of  signal,  noise,  and  signal-to-noise  ratio. 

12? 


C(JK) 

Temperature  values  used  in  labeling  the  tables 
of  signal,  noise,  and  optimum  slit  width. 

ELEMT 

Element  symbol  to  label  a page  of  tables. 

XLIKE 

Wavelength  of  the  line,  in  Angstroms,  to  label 
a page  of  tables. 

XETP 

Description  of  the  flame  type  used  to  label  a 
page  of  tables. 

RATIO 

Description  of  the  fuel  to  oxygen  ratio  used 
to  label  a page  of  tables. 

I 

Index  limit  of  subscript  J identifying  in- 
tensity values. 

M 

Index  limit  of  subscript  K identifying  arbi- 
trary slit  width  values. 

KJ 

Index  limit  of  subscript  JK  which  is  used  to 
select  the  proper  temperature  labels  in 
printing  the  tables  of  signal  and  of  noise. 

L1,L2,L3,IA 

Index  limits  used  to  select  proper  temperature 
values  in  printing  tables  of  signal-to-noise 
ratio. 

N1,N2,N3,N4 

Values  of  temperature  used  in  labeling  the 
signal-to-noise  ratio. 

tables 

XKD 

kd 

XXC 

kc 

XDQ 

Af 

XI D 

id 

XWM 

Wm 

m 

XhT,NL 

Index  limits  controlling  the  number  of  separate 
flame  type  data  decks  to  be  read  into  the  com- 
puter. 

XRD 

Ed 
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XP  \j) 

XIC  ic 

The  composition  of  a complete  set  of  data  decks  for  a 
single  element  and  spectral  line,  assuming  8 values  of  slit 
width  and  the  same  number  of  other  values  as  were  used  in 
the  previous  program,  is  as  follows: 

(Deck  of  values  that  are  independent  of  flame  type) 

I M KJ  ELEMT  XLIXE  SW(1)  SW(2)  SW(3)  SW(h)  SW(5) 

SW(6)  SW(7)  SW(8)  D(l)  D(2)  D(3)  D(4)  D(5)  D(6) 

D(7)  D(8)  D(9)  D(10)  D(ll)  D(12)  D(13)  0(14)  D(15) 


D(16) 

D(17)  0(18) 

D(19) 

D(20) 

D(21) 

D(22) 

0(23) 

D(24) 

D(25)  0(26) 

0(27) 

D(28) 

D(29) 

0(30) 

0(31) 

0(33) 

D(34)  0(35) 

D(36) 

D(57) 

0(38) 

0(39) 

L(40) 

0(41) 

D(42)  D(43) 

d(44) 

D(45) 

D(46) 

0(47) 

0(48) 

D(49) 

D(50)  D(51) 

D(52) 

13(53) 

0(54) 

0(55) 

0(56) 

0(57) 

0(58)  D(59) 

D(60) 

D(61) 

D(62) 

D(63) 

0(64) 

XKD 

XKC 

XDQ 

XI D 

XWM 

XW(1) 

XW(2) 

XW(3) 

XW(4) 

xw(5) 

XW(6) 

xw(7) 

XW(  8) 

0(1) 

0(2) 

0(3) 

0(4) 

0(5) 

0(6) 

0(7) 

0(8) 

XRD 

(control  card  governing  the  number  of  flame  type  decks  to 
be  used) 

XL 

(flame  type  deck) 
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XI  c 

XP 

xi(i) 

XI  (2) 

XI(3) 

XI  (4) 

XI  (5) 

XI  (6) 

XI  (7) 

XI  (8) 

XI(9) 

XI(10) 

XI (11) 

XI  (12) 

XI(13) 

XI (14) 

XI(15) 

XI(16) 

XI(17) 

XI(18) 

XI(19) 

XI (20) 

XI  (21) 

XI(22) 

XI(23) 

XI (24) 

XI(25) 

XI(26) 

XI  (27) 

XI (28) 

XI (29) 

XI (30) 

XI (51) 

xi (32) 

XI  (33) 

XI (34) 

XI (35) 

XI(36) 

XI(37) 

XI(38) 

XI(39) 

XI (40) 

XI (41) 

XI (42) 

XI(43) 

XI (44) 

xi (45) 

XI ( 46 ) 

XI(47) 

XI ( 48 ) 

XI (49) 

XI(50) 

xi(51) 

xi (52) 

XI(55) 

XI(54) 

XI(55) 

XI(56) 

XI (57) 

XI(58) 

XI(59) 

XI (60) 

XI (61) 

XI (62) 

XI(63) 

XI ( 64 ) 

An  example  of  the  form  in  which  the  calculated  values  of 
optimum  slit  width  are  printed  out: 

ATOMIC  EMISSION 
NA  5890  AC/OX  = 2/5 

OPTIMUM  SLIT  WIDTH  IN  CM 

LOG  OP  P TOTAL  -10  -9  -8  

TEMPERATURE 

1800 

2000 

2200 

(Values  of  optimum  slit  width) 


• • 0 
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An.  example  of  the  form  in  which,  the  calculated  values  of 
signal,  noise  and  signal-to-noise  ratio  are  printed  out. 
There  is  room  for  complete  sets  of  values  for  four  tempera- 
ture values  on  each  page. 

ATOMIC  EMISSION 

NA  5890  AC/OX  = 2/5 

SIGNAL  TO  NOISE  RATIO 

SLIT  WIDTH  IN  MICRONS  30.  50.  70 

TEMPERATURE  LOG  P TOTAL 
2400  -10 

- 9 

- 8 

. . (values  of  signal-to-noise  ratio) 


2600  -10 

- 9 

- 8 


The  complete  Fortran  program  is  as  follows: 


DIMENSION  XI (64) ,W0(64 ) ,XW( 8) f SNR ( 64 , 8 ) , C ( 8 J ,SW( 16) VD( 
1128) , S I G ( 64  » 8 ) ,RAC( 64, 8 ) 

1 READ  INPUT  TAPE  5 , 90 , I , M , K J , ELEMT , XL  I NE , ( SW ( K ) , K = 1 , M ) , 
1 { DC  J ) , J=1 v I) 

READ  INPUT  TAPE  5 , 92 , XKD, XKC , XDQ , X ID , XWM , { XW ( K ) , K=i , M ) 
i , (C( JK) , JK=I ,KJ) » XRD 

WRITE  OUTPUT  TAPE  6 , 94 , I , M , K J , ELEMT , XL INE , ( S W ( K ) , K= l , M 

1) »{D(J)»J-1»I) 

WRITE  OUTPUT  TAPE  6 , 96 , XKD , XKC , XDQ , X I D , XWM , ( XW ( K ) , K= 1 , 
IM) , (C( JK) , JK=1,KJ) » XRD 
90  EORMAT  ( 3 1 4 , 2 A5 , 10F 5 . 0/ { 2X , 1 4F 5 . 0 ) ) 

92  FORMAT  (6E12.2) 

94  FORMAT  ( 1H1 , 3 I 4 , 2A5 , 10F 5 .0/ ( 2X , 14F5.0 ) ) 

96  FORMAT  ( 1H0/ ( 2X , 6E 1 2 .2  ) ) 

READ  INPUT  TAPE  5,102, NL 
DO  65  NN= 1 , NL 

READ  INPUT  TAPE  5 , 1 03 , L 1 , L 2 , L3 , L4 , N 1 , N2 , N3 , N4 , XF TP , R A T 
iIO,XIC,XP, (XI ( J) , J= 1, I ) 

WRITE  OUTPUT  TAPE  6 , 104  , L I , L2 , L3 , L4 , N 1 ,N2 , N3 , N4 , XFTP , R 
i AT  1 0 , X I C , XP , (XI (J) , J= 1 , I) 

DO  40  J=l, I 
I KY=  1 
Q=  1 . GE-3 

A=(2.0*XDQ*{  (XKC*XP)**2.0)*(  ( ( ( X I C*XRD ) **2 . 0 ) * ( Q * *4 . 0 ) 
L ) + ( ( ( XWM* ( (XIC»XRD) » * 2 . 0 ) ) + ( X I C*XRD*X I ( J)))*(Q**3.0))) 

2) 

3= ( XKD*  XDG» ( (XKC»Q*( ( X I C»XRD*XWM ) +X I ( J ) ) ) ♦ ( 2 . 0* X I D ) ) ) 
IF (ABSF ( 1.00- A/ 8) -0.05) 39,39,20 
20  I F ( A-B ) 25 , 22 , 22 

22  WRITE  OUTPUT  TAPE  6 , 105 , J , Q , A , 8 

23  GO  TO  39 
25  Q=Q+5.E— 4 

A= ( 2 • 0 *XDQ* ( (XKC*XP)»*2.0)*{ ( ( ( X IC*XRO ) **2 . 0 ) * ( Q# *4. 0 ) 
i ) + ( ( l XWM* ( (XIC*XRD)**2.0) ) + ( X IC*XRD*XI ( J) ))*(Q**3.0))) 
2) 

B= ( XKD*XOQ* ( ( XKC*Q* ( l X I C*XRD*XWM ) +X I ( J ) ) ) + ( 2 . 0* X ID ) ) ) 
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I F ( ABSF ( I. 00- A/ B) -0.05) 39, 39, 27 
27  I KY=  I KY  + 1 

IF  ( IKY-200)20,29,29 

29  WRITE  OUTPUT  TAPE  6 , 106  , J , Q , A , B 

30  CO  TO  39 

39  WO ( J ) =Q 

40  CONTINUE 

WRITE  OUTPUT  TAPE  6 , 10 7 , ELEMT, XL  I NE , XF TP ,RA T 1 0 

DO  50  JK= 1 , K J 

LJ-=J:<*8-7 

JJ=JK*8 

WRITE  OUTPUT  TAPE  6 , 108 ,C ( JK ) , ( WO (J ) , J=L J , J J ) 

50  CONTINUE 
00  60  J = 1,I 
00  60  K= 1 » M 

SIGI J,K)=XKC*XW(K)*XI  ( J ) 

RAC ( J , K ) = { ( XCQ* ( (XKD*(XIQ+(XKC*XW{K)*( { X I C*XRD* ( XW ( K ) + 
1XWM) ) +X I { J )))))+ ( (XKC*XW(K)*{XP*XIC*XRD*(XW(K)+XWM)+(X 
2 P*X  I (J) ) ) )**2) ) ) ** • 5 ) 

SNR( J,K)=SIG( J,K)/RAC( J,K) 

60  CONTINUE 

WRITE  OUTPUT  TAPE  6 , 109 , ELEMT , XL  I NE , XFTP , RAT  10 , ( SW l K ) , 
1 K= 1 » M ) 

LJ=L 1*8-7 
J J=L 1*8 
LL=L2*8-7 
KK=L2*8 
MM=L  3*  8“ 7 


NK-L  3*8 

KM=L4*8-7 

KN’=L4*8 

WRITE  OUTPUT  TAPE  6 , 1 10 ,N1 , { D ( J ) , ( SNR ( J , K ) , K=1 , M ) , J=L J 
1,JJ) 

WRITE  OUTPUT  TAPE  6 , 1 10 , N2 , I D ( J ) , ( SNR ( J , K ) , K=1 , M ) , J=LL 
1 »KK) 

WRITE  OUTPUT  TAPE  6 , 1 10 ,N3 , ( 0 ( J ) , ( SNR ( J , K ) , K=1 , M ) , J=MM 
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1 »NK) 

WRITE  OUTPUT 
1 ,KN) 

WRITE  OUTPUT 
1K=1,M) 

UO  61  JK= 1 ,4 
LJ=JK*8-7 
JJ=JK*8 
WRITE  OUTPUT 
1 = L J » J J ) 

61  CONTINUE 
WRITE  OUTPUT 

1K=1,M) 

DO  62  JK=5 , 8 
LJ=JK*8-7 
JJ=JK*8 
WRITE  OUTPUT 
1=LJ, JJ) 

62  CONTINUE 
WRITE  OUTPUT 

1 K= 1 » M ) 

00  63  JK= 1 , 4 
LJ=JK*8-7 
J J-=  JK*8 
WRITE  OUTPUT 
1=LJ, JJ) 

63  CONTINUE 
WRITE  OUTPUT 

1K= 1 » M ) 

DO  64  JK=5 » 8 
LJ= JK*8-7 
J J= JK*8 
WRITE  OUTPUT 
1=LJ, JJ) 

64  CONTINUE 

65  CONTINUE 


TAPE  6 , 1 10  » N4 , (0(J) » (SNR( J,K) ,K=1,M) ,J  = KM 
TAPE  6, 111 ,ELEMT,XLINE,XFTP,RATIO, (SW(K) , 

TAPE  6, 112, Cl JK), ID(J), ( S IG ( J , K ) , K=  1 ,M ) , J 
TAPE  6, 111 , ELEMT, XL INE.XFTP, RATIO, { SW(K) , 

TAPE  6, 1 12 , C { JK ) , ( 0 ( J ) , (SIG(J,K),K=1,M) , J 
TAPE  6, 11 3, ELEMT, XL INE.XFTP, RAT  10, (SWIK) , 

TAPE  6, 112,C( JK), (D( J ), (RAC( J,K),K  = 1,M)  , J 
TAPE  6, 11 3, ELEMT, XL INE.XFTP, RAT  10, (SW(K) , 

TAPE  6,112,C(JK),(D(J),(RAC(J,K),K  = 1,M)  ,J 
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gg  ru  1 

102  FORMAT  (1112) 

103  FORMAT  (8I6/2A5/ { 6E  12.2  ) ) 

104  FORMAT  ( 1H2 , 8 1 6/4X , 2A5/ { 6E 12 . 2 ) J 

105  FORMAT  ( 1H0 , 4X , 5 1HW-ZERC  OVERSHOT  TRUE  VALUE  BY  LESS  T 
1HAN  5 MICRGNS/4X, 115, 3E12.2) 

106  FORMAT  i 1H0 , 4X, 38HW-ZER0  IS  GREATER  THAN  ONE  MILLIMETE 
1R  /4X,125,3E12.2) 

107  FORMAT  ( 1H1 , 3X , 15HAT0MIC  EM  I SS I CN/4X , 4A5 // /49X , 2 4HGPT I 
1MUM  SLIT  WIDTH  IN  CM/// 5X , 1 4HL0G  OF  P TO TAL , 8X , 3H- 10 , 1 
10X,2H-9, 10X,2h-8, 10X,2H-7, 1 OX, 2H-6 , 10X , 2H-5 , 1 OX , 2H-4 , 1 
30X.2H-3//9X, 11HTEMPERATURE ) 

108  FORMAT  ( 1H0 , 1 IX , F6. 0 , 8E 12 . 2 ) 

109  FORMAT  ( 1H1 ,3X, 15HAT0MIC  EMISSION// 4X,4A5///48X, 2 1HSIG 
INAL  TO  NOISE  R AT I0///4X , 2 1 HSLI T WIDTH  IN  MICRONS, 2X,F 

25.0. 7X,F5.0,7X,F5.0 ,7X, F5.0,7X, F5.0, 7X , F 5 . 0 , 7X , F 5 . 0 , 
37X,F5.0,//4X,24HTEMPERATURE  LOG  P TOTAL) 

110  FORMAT  ( 1H0,7X, I 5/ ( 17X , F5 . 0 , 2X , 8E12 . 2 ) ) 

111  FORMAT  l lhi , 3X, 15HAT0MIC  EM  I SS I ON/ /4X , 4A5/ / /5 1 X , 14HS I G 
INAL  IN  AMPS///4X.21HSLI T WIDTH  IN  MICRONS, 2XF500»7X,F5 
20 , 7X , F5 . 3, 7X , F5 . 0 , 7X , F5 . 0 , 7X , F5 • 0 , 7X , F5 . 0 , 7X , F5 . 0//4X , 
324H TEMPERATURE  LOG  P TOTAL) 

112  FORMAT  { lH0,7X,F6.0/( 1 7X , F 5 - 0 , 2X , 8E 1 2 . 2 ) ) 

113  FORMAT  (1H1.3X, 15HA10MIC  EMISS I0N//4X, 4A5///52X , 13HN0I 
1SE  IN  AMP  ///4X.21HSLIT  WIDTH  IN  MICRONS , 2X , F5. 0 , 7X, F5 

1.0. 7X,F5.0,7X,F5.0, 7X , F 5. 0 , 7X , F5 . 0 , 7X , F5 . 0//4X » 24HTEMP 
1ERATURE  LOG  P TOTAL) 

END 


APPENDIX  4 


Literature  Sources  for  Values  of  Parameters  Used  in  Calcu- 
lation of  Integrated.  Intensity,  Absorption  Coefficient,  and 
Signal-to-Noise  Ratio. 


2 

A area  of  slit,  cm 
a classical  damping  factor 
A^  transition  probability,  sec-1 

B factor  characteristic  of  photo- 
detector surface 

B(P)  partition  function 

E excitation  energy  of  state  u, 
electron  volts 

F focal  length  of  collimator,  cm 
gu  statistical  weight  of  upper  u state 

E height  of  slit 

R^  reciprocal  linear  dispersion,  mp  cm' 

transmission  factor 

minimum  resolving  power  slit  width, 
mp 

factor  to  account  for  broadening 
other  than  Doppler  broadening 


24 

Vickers 

Hinnov1^ 

Corliss  and  Boz- 
man^ 

24 

Vickers 

Landolt-Born- 

stein^ 

Landolt-Bc>rn- 

stein^ 

24 

Vickers 

Landolt-Born- 

stein^ 

24 

Vickers 

24 

Vickers 

24 

Vickers 

Herrmann  and 
Alkemade® 

21 

Poesner 
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